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AN INTEGRAL ANALYSIS OF HEAT TRANSFER DOWNSTREAM OF A 
REARWARD- FACING STEP WITH SMALL COOLANT INJECTION 

By Keiichi Karashima* 

Langley Research Center 

SUMMARY 

This paper presents an approximate analysis of heat transfer downstream of a two- 
dimensional rearward-facing step followed by an infinite flat plate placed in a uniform 
supersonic stream. It is assumed that the thickness of the oncoming boundary layer is 
small compared with the step height and, that the coolant mass is injected into the dead- 
air region through the step face at a small flow rate. The whole flow field is divided into 
three regions: a jet mixing region, a recompression region, and a fully developed turbu- 
lent boundary layer. 

It is shown that if there is no coolant flow, the flow separation accompanying a 
rearward-facing step and the associated recompression provide rather high heating rates 
in these regions when compared with the corresponding flat-plate flow. In the range of 
a very low coolant flow rate, although local heating rate in these high heating regions 
decreases, net heat transfer to the wall in these regions increases slightly because of 
their elongation. 

On the other hand, a large amount of mass transport due to turbulent mixing is 
shown to result in considerable heat -transfer reduction because of growth of the boundary- 
layer thickness in the region downstream of recompression; thus, a favorable compensa- 
tion for the increased heating encountered upstream is indicated. 

It is shown that the overall effectiveness of coolant injection becomes a maximum 
at a very low coolant flow rate, and higher effectiveness can be obtained for smaller 
thicknesses of oncoming boundary layer; thus, an upstream location of the coolant slot 
is preferable. Moreover, examination of effect of coolant slot height reveals that if 
thickness of the oncoming boundary layer is fixed, there exists an optimum coolant slot 
height for which overall effectiveness of film cooling is a maximum. 


iNRC-NASA Resident Research Associate, now at University of Tokyo, Tokyo, 
Japan. 



INTRODUCTION 


Film cooling is used to alleviate severe aerodynamic heating which is encountered 
in many types of engineering problems with high-enthalpy gas flow. It is used, for 
instance, to cool walls of combustion chambers and nozzles in jet engines and rocket 
motors, blades in gas-turbine engines, and nozzle walls of high-enthalpy wind tunnels. 

For entry from earth orbit, film cooling may have the capability of replacing conventional 
ablative heat shields. 

A number of studies have been conducted on the film-cooling problem, mainly by 
experimental methods, over a long period of time; therefore, considerable information 
is available about the correlation between film cooling effectiveness and coolant flow 
rate. (See refs. 1 to 5.) The effect of coolant flow rate on adiabatic wall temperature 
has been established empirically for the case of incompressible boundary- layer flow 
over a heated surface, and analyses for this case are presented in references 6 and 7, 
for example. 

The present approach to the problem is to treat the case of the heat transfer rate 
downstream of a rearward-facing slot where coolant is injected tangential to the surface. 
The wall is maintained at an isothermal condition, for example, by first vaporizing the 
fluid on the inside of the vehicle wall, a condition of interest for some entry vehicles. 

No theoretical approach has been developed for solving this problem. The problem is 
complex, because the concentrated injection at a given location produces an abrupt change 
of boundary layer characteristics, and this in turn causes a complicated interaction 
between the internal dissipative flow and the external isentropic supersonic flow. More- 
over, if the internal dissipative flow is turbulent, even theoretical knowledge of the effect 
of compressibility on the flow is meager. 

The purpose of the present paper is to demonstrate the importance of turbulent 
mixing in determining the characteristics of dissipative flow downstream of the point of 
coolant mass injection and to formulate this concept in quantitative terms to clarify the 
effect of various aerodynamic parameters. The resulting equations are solved numeri- 
cally for the isothermal wall condition, and the solutions are examined in terms of the 
gross heat transfer characteristics downstream of the slot. 

To formulate and solve the problem, several assumptions are introduced which 
limit the applicability of the results. First, based on the experimental results of ref- 
erence 4, the problem is related to that associated with the base pressure behind a 
two-dimensional rearward-facing step through which mass is injected into the dead-air 
region, so that the configuration analyzed is a rearward facing two-dimensional step 
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followed by an infinite flat plate. Further, the oncoming boundary layer thickness is 
assumed to be small compared with the step height. Also, the coolant flow rate is 
assumed to be small, so that the velocity in the dead-air region can be neglected. Other 
assumptions, which do not reduce the scope of the problem being solved, are introduced 
in the theoretical development as required. 

SYMBOLS 

a speed of sound 

C length of plate lip 

Cf skin-friction coefficient 

Ch,oo energy transfer coefficient defined by equation (87) 

C n constant of integration 

Cq heat-transfer coefficient 

Cj constant of proportionality 

Cp specific heat at constant pressure 

D height of step 

e internal energy 

F,G functions defined by equations (16b) and (16c) 

Fj,F 2 functions of k defined by equations (68) 

f function 

H stagnation enthalpy 

h static enthalpy 
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momentum flux defined by equation (47) 


J energy flux 

K mixing coefficient 

k D> k M constants of proportionality with respect to boundary -layer thickness 
L length of separated region 

l width of recompression region 

M Mach number 

]VI coolant flow rate defined by equation (33) 

m reduced mass flux defined by equations (58) 

m mass flux defined by equation (46) 

in nondimensional mass flux defined by equations (74) 

Npr,t turbulent Prandtl number 

N Re Reynolds number, n^C/v^ 

N Re d Reynolds number, UoqD/^oo 

Nr 6> l Reynolds number, u-^L jv^ 

N R e oo Reynolds number, u^x/v^ 

P function of k and We defined by equation (62) 

p pressure 

Q heat input defined by equations (35) 

q heat -transfer rate ' 
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q 

q w 

R 

s 

T 

T 

u,v 

u 

V 

W p 

W s 

w 

X 

x,y 

X 

ym 

z 

a >fi 

y 


mean heat -transfer rate 

heat-transfer rate on a flat plate (Van Driest' s solution) 
gas constant 

surface of control volume 
temperature 

mean temperature across viscous layer in recompression region 

velocity components 

mean velocity defined by equation (49) 

resultant velocity 

work done on gas by pressure, defined by equations (35) 
work done on gas by shear, defined by equations (35) 
reduced velocity defined by equations (58) 
constant of proportionality pertinent to scale factor 
coordinates 

nondimensional distance defined by equations (74) 
mean thickness of mixing layer 

p 

parameter defined by — 



constants defined by equations (90) 
ratio of specific heats 
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thickness of boundary layer 


displacement thickness 

nondimensional thickness of boundary layer defined by equations (74) 
angle between external flow direction and X-axis in jet mixing region 
similarity variable defined by equation (8) 

local effectiveness of coolant injection in region of fully developed 
turbulent boundary layer 

overall effectiveness of coolant injection 

flow deflection angle of external stream 

momentum thickness 

parameter defined by equation (50) 

wall enthalpy ratio defined by 

viscosity 

kinematic viscosity 

coordinates in corresponding incompressible flow 
density 

mean density across viscous layer in recompression region 
scale factor 
shearing stress 


parameter defined by equation (64) 



r 


x/s stream function 

co power index for viscosity-temperature law 

Subscripts: 

0 stagnation condition of oncoming flow or conditions without coolant flow 

1 external flow conditions in jet mixing region 

2 flow conditions at point of injection 

b conditions at base of step without coolant flow 

c conditions of coolant flow at slot exit 

d conditions on dividing stream line 

e external stream conditions in recompression region 

i conditions in corresponding incompressible flow or conditions at exit of 

jet mixing region 

j conditions in jet mixing region 

m conditions in mixing region 

p flat plate 

R conditions in recompression region 

r conditions at entrance of recompression region 

s stagnation conditions 

u conditions at upper jet boundary 

w conditions at wall 
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oo 


free stream approaching the step 


Superscript: 
n integer 

A bar over a symbol denotes a nondimensional value; a tilde denotes a mean value. 
An asterisk denotes that the value has been normalized by the reference values. Primes 
denote differentiation with respect to £. 

PHYSICAL MODEL 


As a flow model to be analyzed, consider a two-dimensional rearward-facing slot 
followed with an infinite flat plate placed in a uniform supersonic stream. To simplify 
the analysis, the coolant flow is assumed to consist of the same species as the main 
stream and to be injected at mean subsonic velocity u c and density p c at the slot exit. 

A schematic illustration of the characteristic flow pattern is shown in figure 1. 
Since the oncoming flow is supersonic, a Prandtl-Meyer fan emanates from the plate lip 
to deflect streamlines inward, and a separated boundary layer spreads itself downstream 
like a free jet. Inasmuch as the flow must become parallel to the flat plate far down- 
stream, a recompression shock wave will make an appearance at some distance down- 
stream of the slot, to deflect streamlines outward, and the viscous layer again attaches 
to the plate. When this simplified flow model is specified, the whole flow field can be 
divided into three regions. The first is a jet mixing region behind the point of 


Oncoming boundary layer 

Prandtl-Meyer fan 



K Jet mixing - 


i Fully developed - 

f turbulent 

I 

^ Recompression 


Figure 1.- Illustration of overall flow pattern. 
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separation, in which mixing between the external isentropic stream and internal dissipa- 
tive flow determines the flow characteristics. The second is a recompression region 
next to the jet mixing region, in which an interaction between a shock wave and the viscous 
layer is the essential physical feature. Finally, there is an attached boundary layer along 
the flat plate downstream of the recompression region. 

It is of primary interest to specify whether the separated viscous flow in the jet 
mixing region is laminar or turbulent. From the experimental point of view, it has 
already been reported in references 8 and 9 that the flow separation associated with an 
abrupt change of wall contour such as a rearward-facing step tends to induce transition a 
short distance downstream of the separation point even if the oncoming boundary layer is 
laminar. Of course, if the Reynolds number is sufficiently small, the separated viscous 
layer may be laminar. However, for practical applications, it seems realistic to assume 
that the separated viscous layer in the jet mixing region is turbulent except in the vicinity 
of the separation point, even with a laminar approaching boundary layer. Consequently, 
the viscous layer downstream is expected to remain turbulent throughout, although there 
is neither a rigorous theoretical proof nor an empirical criterion for this statement. If 
the oncoming boundary layer is turbulent, there is no doubt about this assumption. 

TURBULENT JET MIXING REGION 

The essential feature of the flow in this region is viscous mixing following flow 
separation at the edge of the plate lip. Therefore, mass, momentum, and energy trans- 
port through mixing between the external isentropic stream and internal dissipative flow 
play an important role in determining the flow field. This statement means that in con- 
trast to the usual Prandtl boundary-layer theory, the external flow cannot be regarded as 
a known datum for calculating internal dissipative flow; rather, the development of dis- 
sipative flow itself helps to determine the external flow. In fact, the experimental results 
presented in reference 4 show that the strength of the Prandtl-Meyer fan decreases 
fairly rapidly with increase of the coolant flow rate. Consequently, the Mach number of 
the external flow changes with the coolant flow rate. This condition, in turn, changes the 
strength of the associated recompression shock wave which is again felt by the jet mixing 
layer through the dead-air region. Thus, the injection of coolant flow through a rearward- 
facing step-like slot involves a complicated flow mechanism of interaction in the jet 
mixing region. 

Analogous to the usual base pressure problem, the pressure behind a coolant slot 
may be expressed as a function of oncoming flow Mach number, thickness of oncoming 
boundary layer at the separation point, height of coolant slot and coolant mass flow rate, 
and so forth, such as 
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( 1 ) 


( P c u c \ 

p i ” f \ Moo,5oo,D ’ pZ u^» * * 7 

Although much effort has been devoted to this problem, there does not exist any general 
analytical approach suitable for estimating the pressure rigorously. The main reason 
is the difficulty of treating the viscous interaction mentioned in the previous paragraph, 
since, both in theory and in experiment, information on detailed flow structure in the 
dead-air region is meager. Besides, if the flow in the jet mixing region is specified as 
turbulent, another difficulty is that no fundamental equation is available. This problem, 
in turn, makes it necessary to introduce some empirical information or assumptions in 
any theoretical approach to the jet mixing region in general. Thus, it seems hopeless to 
develop an exact approach to the problem under consideration. 

In spite of these difficulties, however, there is some empirical information pertinent 
to flow phenomena in the jet mixing region. Chapman et al. (ref. 9) pointed out from their 
experiments that an isentropic supersonic flow outside the dissipative region associated 
with flow separation behind a rearward-facing step can be considered to be approximately 
uniform. Exactly speaking, there must be some pressure gradient in the mixing region. 
However, intuitively, this pressure gradient is very small compared with that in the 
recompression region; therefore, an assumption of constant pressure at the edge of 
mixing layer is reasonable. 

Moreover, if the thickness of the mixing layer is small, the usual boundary- layer 
approximations may be applicable to the jet mixing region. Based on this approximation, 
reference 10 presents a momentum equation for incompressible turbulent mixing by use 
of Prandtl's mixing length hypothesis and obtains a similar solution for the velocity dis- 
tribution across the mixing layer. This solution, if the value of a scale factor is chosen 
properly, is shown to be in good agreement with experimental results presented in ref- 
erence 11. Furthermore, the experimental results of reference 8 for compressible tur- 
bulent jet mixing indicate that the theoretical velocity profile presented in reference 12 
agrees well with the experiment, but the value of the scale factor is seriously influenced 
by the Mach number of the external flow. In reference 8, it is further commented, in 
spite of the similar solution being derived on the basis of a concept of zero thickness of 
the oncoming boundary layer, that the similar velocity profile agrees well with experi- 
mental data, even in the case of nonzero thickness of the oncoming boundary layer. 

It must be noted that although the experimental evidence may justify the applicability 
of the assumption of a similar flow field even in the case of nonzero thickness of the 
oncoming boundary layer, this assumption does not imply that the oncoming boundary 
layer has no influence on the flow characteristics in the jet mixing region. In fact, as 
will be shown later, the basic equations governing this region involve the Mach number 
of the external flow explicitly. Since the Mach number and, consequently, the associated 
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pressure in the jet mixing region are influenced by the thickness of the oncoming bound- 
ary layer and, in particular, this trend is strongest for the laminar approaching boundary 
layer, the flow characteristics in the jet mixing region depend clearly upon the oncoming 
boundary-layer thickness. 

This empirical evidence suggests the possibility of developing an approximate 
analytical approach to compressible turbulent jet mixing by using a concept of similar 
solutions, even in the case of nonzero thickness of the oncoming boundary layer. How- 
ever, it must be noted that from a mathematical viewpoint, this similar solution is an 
asymptotic form of an exact solution far downstream of the separation point and, there- 
fore, it is not expected that it can predict the flow field accurately just downstream of 
the separation point. 

The following section brings out the importance of the transport of mass, momentum, 
and energy from the outer stream to internal dissipative flow in determining the flow pat- 
tern in the jet mixing region and formulates this concept in quantitative terms on the 
basis of a simplified flow model. The basic assumptions are summarized as follows: 

(1) Oncoming boundary layer is either laminar or turbulent and is thin compared 
with the step height, whereas the flow in the jet mixing region is turbulent except in the 
vicinity of the separation point where the oncoming boundary layer is laminar 

(2) The pressure gradient in the jet mixing region is negligible 

(3) The turbulent Prandtl number is constant 

Analysis 

Before proceeding with the analytical development, a coordinate system related to 
the mixing region and a condition inside the dead-air region must be specified. Unfortu- 
nately, there cannot be found, a priori, any definite intrinsic boundary within the mixing 
layer to which a physical coordinate system is referred. However, if the assumption of 
a similar flow field is introduced, it provides a characteristic streamline in the mixing 
layer along which physical properties are kept constant. Therefore, in the present 
approach, the x-axis of the coordinate system (x,y) is taken along this characteristic 
streamline with its origin at the point of separation, the y-axis being normal to the 
x-axis. However, it must be noted that the location of the characteristic streamline 
relative to the wall is not known a priori. 

A condition inside the dead-air region is another point of interest. If coolant mass 
is injected into the dead-air region through a rearward-facing slot, the coolant flow must 
have, on the average, a finite nonzero velocity at the slot exit. As it flows downstream, 
mixing occurs between separated boundary- layer flow and the coolant flow by the 
"scavenging" effect. However, since the wall exists, growth of the mixing layer is 
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restricted and the coolant flow near the wall will be decelerated by viscous forces far 
downstream of the injection point. Thus, when the end of the jet-mixing region is 
approached, it can be expected that the mixing has almost finished and the separated 
boundary layer is ready for reattachment. 

Therefore, near the end of the jet mixing region, the mixing layer has grown to the 
extent of almost the entire cross section of the dissipative flow, and the velocity profile 
must be nearly similar to that at the separation point. 

On the other hand, since detailed flow structure in the dead-air region is ignored, 
it is intended in the present approach to estimate overall aerodynamic characteristics 
pertinent to the jet mixing region by considering a mass, momentum, and energy balance 
within a control volume bounded by the cross section at the entrance and exit of this 
region. For this purpose, it is more important to evaluate velocity and enthalpy profiles 
across the dissipative layer in a realistic form at the exit of the jet mixing region rather 
than at the intermediate points. For this reason, u = 0 is assumed as a boundary con- 
dition inside the dead-air region to obtain a velocity profile of the separation type. 

By Prandtl's mixing length hypothesis, local shearing stress and heat-transfer rate 
may be expressed as 

r = e(x) p 


q 


e(x) 3h 
Npr,t 


where e(x) is the eddy kinematic viscosity which is defined as 


e(x) = CjUjX 

and where c^ indicates a constant of proportionality. Thus, the differential equations 
for mass, momentum, and energy conservation may be written within the accuracy of the 
boundary-layer approximation as 


9pu t dpv _ Q 
dx ay 


9u Oil /x d f 

pu & + pv l>r e(x) 5y(' > 


9u\ 

zy) 


9h 


ah 


pu ^ +pv ^ = 


e(x) 9 
Npr,t ®y 



+ £<x),p t 


where all physical quantities indicate time-mean values. 


(2) 

(3) 

(4) 
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The continuity equation may be accounted for by introducing a stream function 




dxly 

_ r* s pu 

Introducing the nondimensional expressions 


x * 

L 


* u 
u* = — 
u 


h* = h 


1 

L_ 

hi 


P* =£- 

P Pi 

V = 


* 


N 


UjL 

UjL 


Re,L 

and rewriting equations (3) and (4) then yields 


where 


8u* _ 

e*x* -4 


8u ^ 

\ 

ax* 

a \p 

\h* 2 

at//* 

) 

ah* _ 

X* 

8 i 

^u* 

ah’ 1 

ax* 

N Pr,t 

at//*' 

\h* 2 

at// 

e* = c 

l N Re,L 





+ e*x* u 


* /^.*\ 2 u ^2 


*2 


f 9u 

3t//* 


(5) 


(6) 


(7) 


and where a relation p* = 1/h* has been employed. 

To obtain similar solutions to equations (5) and (6), a similarity parameter based 
on the streamline coordinate is introduced as an independent variable and is given by 

? = a l ( 8 ) 


where denotes a scale factor. Thus, the momentum and energy equations may be 

reduced to 
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( 9 ) 


« _ du* _ d / u* du*\ 
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d /u* dh 


N 


Pr,t d ?lh* 2 d ? 


+ (r - 1)M 1 


2 u* /du*f 

i^2\ d ?y 


( 10 ) 


a, 2 — 

1 c,N 


1 Re,L 


The boundary conditions are given by 

u* = h* = 1 
u* = 0 

h* = *4 



If the wall enthalpy is denoted by H w , then, 


h 


* 

w 


Hw _ H W / y - 1 
hj " Hoo\ + 2 


M, 


(ID 


( 12 ) 


(13) 


In the case where the turbulent Prandtl number equals unity, the energy equation 
can be readily integrated analytically with the assumption that local static enthalpy is a 
function of u* only. The result is (Np r ^ = 1.0) 


h* 


y - i 

2 


Mj 2 u * 2 


+ 



(14) 


It remains to discuss the value of the scale factor which must be determined 
empirically. In the case of incompressible jet mixing, it is shown that if a similarity 
parameter £ is defined by the use of an intrinsic coordinate system x,y such as 

then good correlation between the theoretical velocity profile presented in reference 10 
and the experimental data of reference 11 can be obtained for a - 13.5 (ref. 12). In the 
case of compressible jet mixing, it is further shown (ref. 8) that the experimental velocity 
profile across the mixing layer at = 1.6 can be well correlated with the theoretical 
profile proposed in reference 13 for incompressible jet mixing, if the value of a is 
chosen to be 18. On the other hand, a formula for o for compressible jet mixing which 
accounts for Mach number effects can be derived by fitting an equation to the data of ref- 
erence 11 at = 0 and reference 8 at = 1.6. This scale factor is given by the 
equation 
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a= 13.5 + 


On the other hand, since the similarity parameter £ used in the present approach 
is based on a streamline coordinate, the associated scale factor becomes different 
from that for the intrinsic coordinate system. However, it can be shown from refer- 
erence 8 that the constant of proportionality involved in the eddy kinematic viscosity is 
inversely proportional to a 2; that is, 



From this relation together with equation (11) it is found that 
CT W N Re,L = -j= “ CT 


Therefore, the scale factor is given by the equation 


a l\/ N Re,L = x (* 3 - 5 + 6 ' 16M 1 < 15 ) 

where X denotes a constant of proportionality. Unfortunately, the value of X is not 
known rigorously. As will be shown in the subsequent section, the value of X depends 
strongly upon the value of a mixing coefficient which seems to play the most important 
role in determining the aerodynamic characteristics in the jet mixing region. Since there 
exists neither theoretical nor empirical information on rigorous evaluation of the mixing 
coefficient, it is assumed for the present that X is unity. To demonstrate the effect 
of X on flow characteristics, an example of the heat-transfer rate is calculated later 
for X = 2 and the results are compared with those for X = 1 . 

Equations (9) and (10) can be integrated by using an iteration scheme. Let 

W ( 16 a) 

h*^ 

^ff =F(?> (16b ' 

jfrfF- 0 ® (i6c » 

Then solutions for u* and h* satisfying the given boundary conditions may be obtained, 
respectively, as (see appendix) 


u 


* 




(16d) 
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F Npr ,t 


i: 


d? 


d? 


-2 p p? F N Pr,t ( -4 F ^" N Pr,t * 

- ^ N Pr,t^ " D M 1 J J y d ?7 d£ + 1% 


( 17 ) 


where 


F(0 = exp 




Numerical calculation is carried out by evaluating f(£), first, from assumed pro- 
files for velocity and static enthalpy as a first-order approximation. Then, by use of the 
f(£) value thus obtained, the first-order velocity and enthalpy profiles are calculated 
from equations (16) and (17), respectively. With these values, the second-order approxi- 
mation of f(£) is obtained and the same calculation procedure is repeated until solutions 
for u* and h* converge sufficiently. 

Mass and Momentum Balance at Exit of Jet Mixing Region 

Although the rate of spread of the mixing layer itself is known from the value of the 
scale factor, the location of the characteristic streamline (£ = 0) relative to the wall is 
not known a priori. However, it may be determined by considering the transport of mass 
and momentum from the external stream to the internal dissipative flow. 

It is clear that in the jet mixing region, there are two physically important bound- 
aries other than the characteristic streamline. One is the upper jet boundary which 
separates the region of internal dissipative flow from the external isentropic supersonic 
stream. Another is the so-called "dividing streamline" which separates the mass of 
fluid originally in the oncoming flow from the coolant flow. However, at a free boundary 
such as jet boundary, any physical quantity in the dissipative region must approach 
asymptotically the value of the outer flow; thus, there is no intrinsically distinct boundary 
between the dissipative region and the external flow. Because of this, the upper jet bound- 
ary is defined in the present approach as the loci of points where u* = 0.99. 
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Simple sketches of the characteristic flow pattern shown in figure 2 may be helpful 
for subsequent development. From conservation of mass it is easily shown that the mass 
of fluid entering across segment AB (fig. 2) must be equal to that leaving across segment 
CF at the exit of mixing region (x = L or x* = 1); that is, 

(J^ u pu dyj - Pi u i(y u - - 5 i) + Pi u i( a i - V) 

The first term on the right-hand side of the preceding equation indicates the mass of 
fluid entrained in the control volume ABCF from outer inviscid flow through mixing, and 



e 

— f y* 


m 


«^3 


Streamline direction 
of external flov 


" xU = 0) 


Vn U \ \ \ \ \ \ \ \ >A \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ V 

Dividing streamline 


\ uW\W\n\\ 
X = L 


(a) 



(b) 


Figure 2.- Characteristic flow field near exit of jet mixing region. 


17 



the second term denotes the mass of fluid originally existing in the oncoming boundary 
layer. Thus, the preceding equation may be reduced by use of nondimensional expres- 
sions to 

1 ?U $ « - + "i tfw TT (18> 

From the definition of the dividing streamline, the mass of coolant flow entrained in the 
mixing region must exist in a region between the dividing streamline and the wall (fig. 2); 
that is, 



which can be reduced to 


“u * 5. 


Vji 

a W N Re,L 


( 19 ) 


On the other hand, the momentum of fluid involved in the viscous layer at x = L 
may be expressed as the sum of the momentum entrained in the viscous region from 
external flow by mixing, and the work done on the gas due to shearing stress at the wall; 
that is, 


f U pu 2 dy = P 1 u 1 2 (6 1 - 6j* - 6^*)+^ 

C Yu pu dy - p^fsj - 

J y i x 

_ ^d 


+ u c J Yd pu dy - r W)j L 


where t ,, 7 ^ denotes mean shearing stress on the wall. The first term on the right-hand 
side of the preceding equation indicates the contribution from the oncoming boundary layer 
and the second and third terms denote momentum entrained from external isentropic flow 
and coolant flow, respectively. The last term is the work done by gas on the wall. By 
use of nondimensional expressions, the equation may be reduced to 


■Cu 

5i 


u*d? = c u - u c *?i - (l - u c % - 


0 i_ 

L 


+ 2 


(20) 


where 


C f,j = T 


JLA. 


2 p l U l 
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Moreover, from the geometry, one can obtain an approximate expression such as 
(see fig. 2) 

D * L sin(0 + e) + | y 4 1 « L (0 + e) +| y A | 

where 


tan e 



1 

^^RejL 




yi 


L 

CT l(/ N Re,L 




and hence 


where 




p? h * 

J(?) = TT d? 

J 0 u 

Thus, elimination of and L from equations (18), (20), and (21) gives 

U(Sm) 
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(22) 
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u c*(? u " k) 


h n 6 l 0 1 
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CT lV N Re,L 0 " J (?i)] + \ CT W N Re,L Cfj 


(23) 


In the case of no coolant flow, where u c * = 0 and = C^, equation (20) is further 
reduced with the aid of the mass balance equation (eq. (18)) to 


»* - ?u - ?d - -,r[ c 'i'l ii R^L e + J ( ?m ) • J ( 5 <0] • \ J 

and the mass balance equation may be reexpressed by use of equation (22) as 

?d - ?u - J(?u) + (l + %)■ J(Cm) + - J(Ci)] 


(24) 


(25) 
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In the preceding equations, 6j* and 0j* indicate displacement thickness and 
momentum thickness of the entrained boundary layer just downstream of Prandtl-Meyer 
fan and can be connected with those of the oncoming boundary layer when mass and 
momentum conservation of boundary-layer flow through the Prandtl-Meyer expansion is 
assumed; that is, 

Pl U l( 6 l - 6 1*) = Poo U oo( 5 oo - Sj) 

Pl U l 2 ( 6 l “ 6 1* " 0 1*) = P=o U oo 2 ( 5 oo - 6 =o* - e J) 

Since the state of the boundary layer will not change through a Prandtl-Meyer fan irre- 
spective of whether the oncoming boundary layer is laminar or turbulent, 6j* and 
j 6^ will not change much from S^*/ 6 m and 0^/fi^, respectively. Thus, rather 
simple conservation equations such as 


Pi Vi* = Poo U oo 5 oc* 

Pl U l 2f9 i* = Poo u oo 20 co* 


can be applied. Therefore, 


6 1 PooUoo 5 0 




D P i u i D Mi 
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01* p oo u oo Soo* M c 
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2 + (y - l)Mj J 
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(26) 


Now 5oo* and 6^* can be expressed in terms of Reynolds number referred to the 
oncoming flow conditions for either laminar or turbulent approaching boundary layers by 
equations such as 


5 * c 

oo \r ^ r-7 

D " d 
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(Equations continued on next page) 
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Equations (24) and (25) can be considered as simultaneous equations with respect to Z 
and where n = 2 for the laminar oncoming boundary layer and n = 5 for the 

turbulent oncoming boundary layer. In the preceding equation K D and K M are 
constants of proportionality associated with displacement thickness and momentum thick- 
ness, respectively. Their values are, of course, changed depending upon whether the 
oncoming boundary layer is laminar or turbulent. 

Thus, in the case of no coolant flow (u c * = o), Z is known for given conditions if 
j^(= £j) is known. Since the oncoming boundary layer is not influenced by coolant mass 
injection, the value of Z must be kept constant for the coolant flow case. Therefore, 
and are known for each coolant flow case and with these values J(? m ) 

and are determined by using equations (23) and (25) if is known. The procedure 
for this calculation is described in detail in a later section. In any case, it seems neces- 
sary, first of all, to know the value of t^. Unfortunately, solutions appropriate to the jet 
mixing region do not give any information on therefore, it must be determined from 
recompression conditions. As is well-known, since a complicated interaction between a 
shock wave and a turbulent boundary layer is the essential physical feature in the recom- 
pression region, it seems hopeless to develop any rigorous approach to this interaction 
phenomenon. To avoid this difficulty, reference 14 proposes an intuitive criterion for 
penetration of fluid particles through the recompression region. This criterion may be 
expressed in the statement that at the entrance of recompression region, the fluid must, 
at least, have a level of mechanical energy equal to the static pressure behind the shock; 
as a result, the recompression to that static pressure is possible by complete conversion 
of the kinetic energy isentropically. Although there is no rigorous proof for this state- 
ment, this criterion could be used. However, this method underestimates the value of Ci, 
since the effect of viscous dissipation in the recompression region is neglected. Because 
of this underestimation, in the present approach will be evaluated from a condition of 
mass flow conservation between the jet mixing region and the recompression region, and 
will be discussed briefly in a later section. For this purpose, it may be helpful to present 
a mass-flow relation at the exit of the jet mixing region, 



(28) 


21 



The fact that equations (23) and (25) involve D and 0 j*/d clearly indicates 

that .physical properties in the jet mixing region, such as pressure and the length of the 
separated region, depend upon the ratio of oncoming boundary-layer thickness to coolant 
slot height. This dependence, in turn, means that such properties can be correlated to 
the parameter Z as in the case of the usual base pressure problem, and thus the cor- 
relation found in reference 15 is confirmed. 

The value of mean skin-friction coefficient Cf j will now be considered. Since 
detailed flow structure near the wall is ignored, Cf j cannot be estimated in the con- 
ventional way; therefore, in the present approach, it is evaluated by use of an approximate 
semiempirical method. For this purpose it is convenient to introduce a concept of a 
mixing coefficient originally proposed in reference 16 which is defined as 

dm- 

ir =K i p i u i 


where Kj is the mixing coefficient and mj denotes mass of fluid entrained in the 
mixing layer from external isentropic flow. The left-hand side of the preceding equation 
may be expressed by 


dm i d f y u ^ p l u l ? u 

a = t- j pu d y = — f — 

Therefore, 

? 

Ki = u (29) 

°i\f^Re^L 

If the parameter <fi is defined for a general state of the viscous layer as 


♦ = . _a . 

2K(1 - k) 

where Cf and K denote local skin-friction coefficient and mixing coefficient, respec- 
tively, and k is a parameter defined by 


5 - 6* - 0* 
6-5* 


then, Cf j may be given in a form 

2 ? (1 - Ki\ „ 

t,i - ^K 1 - *i)*j = - 1 *i 

a iV N Re,L 


C f. 


(30) 
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where Kj and 0j denote their mean values in the jet mixing region. In the case of an 
incompressible turbulent boundary layer on a flat plate, it is easily shown that 0 = 1. 

In general, for any other boundary-layer state, there exists no useful method to esti- 
mate <p. The same difficulty is found in the evaluation of Xj. For this reason, /Cj is 
assumed to be replaceable by the value obtained at the entrance of the recompression 
region; that is, 



On the other hand, the value of seems to depend mainly upon local dynamic 
pressure and wall enthalpy. However, since it seems almost impossible to obtain a func- 
tional dependency of 0j on dynamic pressure and wall enthalpy analytically, the base 
pressure behind a two-dimensional rearward-facing step (u c * = 0) was calculated for 
= 1.5 and 2.0 under the assumption that the base pressure is rather insensitive to 
wall enthalpy. Several sets of parameters ^Moo,Mj,0^, which can be correlated to 
the experimental data shown in reference 15, were found by use of a trial-and-error 
method. The values of these parameters thus obtained correlate through an equation 


where 



S(M) = | yM 2 



r - i 
2 


z 



(31) 


Thus, the final term on the right-hand side of equation (23) or equation (24) is given as 




(32) 


It must be noted that equation (31) is an approximate empirical formula for any future 
rigorous evaluation of the mean-skin friction coefficient in the jet mixing region will 
improve the present approach. 


From the foregoing analysis, physical properties in the jet mixing region such as, 
for example, the pressure, can be expressed as a function of oncoming flow Mach num- 
ber, ratio of oncoming boundary-layer thickness to slot height, and coolant flow rate; 
that is 


p i 




p c u c 
’ Pco U °o’ 



and this equation is an alternate form of equation (1). 
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However, it should be noted that the present approach involves the value of £ u 
explicitly and obviously depends upon the definition of the upper jet boundary. This con- 
dition clearly indicates that the solution to the jet mixing region depends, to a certain 
extent, upon the definition of £ u . Unfortunately, there is no rigorous way to avoid this 
uncertainty. The present approach could be improved with a rigorous definition of £ u 
as well as Cf j. 

It is convenient for subsequent development to introduce a coolant flow rate, nor- 
malized by the oncoming flow conditions, which is defined as 


This equation can be reexpressed by use of local flow conditions in the jet mixing region 
in a form 


M = p c *u c * 


Moo 

2 + (y - 1)M 1 2 

Mj 

_2 + (y - 


y+ 1 

2(y-l) 


(33) 


Energy Balance and Heat-Transfer Rate 

Estimation of heat -transfer rate in the jet mixing region is another important prob- 
lem. As has already been mentioned, the detailed flow structure near the wall is ignored 
in this region; therefore, conventional methods cannot be applied to calculate the local 
heat -transfer rate to the wall. However, it is possible to evaluate a mean heat -transfer 
rate by considering an energy balance within a control volume in the mixing region. 

Consider a control volume enclosed by a surface S which consists of the charac- 
teristic streamline AB, the entrance DA, and the exit BC of the mixing region, and the 
wall CD. (See fig. 3.) Since flow is assumed to be steady, the sum of heat input Q 
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Figure 3*- Control surface for energy balance. 
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across S, the increase in internal energy and kinetic energy within S, the work W p 
done on the gas through pressure normal to S, and the work Ws done on the gas through 
shearing stress on S must be zero; that is, 


Q + E + W p + W s = 0 
Each term is given as follows: 


dS 


E = II( e + f) pVn ds = II( h + t) pV « ds + II pV * ds 

Wp = jj P(-Vn) dS 
S 

- If TV * 


W s = \ \ tv* ds 
' S 


(34) 


(35) 


where the subscripts n and t denote the outward directed normal component and the 
tangential component to S, respectively. Before the estimation of each integral term, 
it is helpful to clarify conditions and assumptions imposed on S, which may be summa- 
rized as follows: 

(1) V n = 0 along AB and CD 

(2) Vt = 0 along CD 

(3) Vt is negligible along BC and DA and 

(4) is negligible along BC and DA 
on 

With these conditions the first integral can be expressed as 



dx - 



where q w j denotes mean heat -transfer rate to the wall. This equation is further 
reduced to 


. P l u l h l L dh«\ , . L 
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where the first term on the right-hand side of the equation indicates heat input into the 
control volume across AB due to conduction. The sum of the second and the third 
integrals in equations (35) may be expressed as 

pO pO 

E + W P = P c u c H c D + J (puH) x _ 0 dy - j (puH) x=L dy 

yd y i 

The first and the second terms on the right-hand side of the preceding equation indicate 
energy inputs into the control volume due to coolant flow and oncoming boundary -layer 
flow, respectively, and the last term denotes the energy going out of the control volume 
across BC. The last integral in equations (35), which indicates heat generation due to 
viscous effects, may be written in a form 


Wc 


- p r«W P. a ax . (r - 

Jo L ^=0 crj^N^; \h* 2 « 


,*2 


/?= 0 


Therefore, equation (34) can be expressed as 
‘ L [ ~ e(x) 

N Pr ,t " 9 y 
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Jy=0 


f0 

dx - q w ^L + p c u c H c D + j (puH) x=Q dy 

yd 


- I y ° (puH) X =L d ? + J 0 L [ £ < x > pu §] y=0 = 0 


(36) 


On the other hand, mass continuity within the control volume gives a relation 


P c u c D 




dy = 


o 


(37) 


where the first and the second terms on the left-hand side of equation (37) indicate mass 
fluxes entrained in the control volume from coolant flow and oncoming boundary-layer 
flow, respectively, and the last term denotes mass flux going out of the control volume 
and across BC. Multiplying equation (37) by Hoo and subtracting the result from equa- 
tion (36), (H) x _q = Heo being assumed and the nondimensional expressions being intro- 
duced, then yields the mean heat -transfer rate in the jet mixing region as 


jgii 

P lUjhj 


L^Re.L^ 




2/u* 2 du*^ 


N Pr,t d? /?=0 
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(38) 
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In the case of Np r t = 1.0, this equation can be further simplified with the aid of 
equation (14) to 



H 
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OO 


OO 


-H c 

-H w 


( ? d - 5i) 


(39) 


where the final term in the bracket of equation (39) indicates that the coolant mass reduces 
the heat -transfer rate as a heat sink. 


If the mean heat -transfer coefficient C Q . 

VL 


in the jet mixing region is defined as 



p^UooHoo 


(40) 


then 

C = P 1 U 1 qw >3 
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(41) 


In contrast to conventional boundary-layer transpiration cooling, in which the Mach 
number of the external flow is rather insensitive to coolant flow rate, with coolant injec- 
tion the external flow Mach number in the jet mixing region is fairly sensitive to the 
coolant mass injection. Furthermore, this characteristic is an unfavorable factor in the 
reduction of heating rate in the jet mixing region, since decreases with an increase 

in the coolant flow rate. This decrease, in turn, tends to increase Cq j, as is seen in 
equation (41), and this fact may be largely responsible for the marked difference in heat- 
transfer rate between separated flow and attached boundary-layer flow. 

Despite the unfavorable factor just mentioned, the film cooling of the present type 
is still effective in reducing heat transfer in the region farther downstream of the sepa- 
ration point because a large amount of mass transport, induced by the strong action of 
jet mixing, can certainly be expected from the external flow to the internal dissipative 
flow. This process increases the thickness of the boundary layer downstream and reduces 
the heating rate to some extent. 

RE COMPRESSION REGION 

One of the most interesting phenomena in fluid mechanics, which is, in an exact 
sense, still unsolved, is the interaction between a shock wave and a viscous layer. This 
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interaction can be observed in many practical problems of supersonic flow, such as base 
pressure behind blunt-based bodies and the flow separation associated with recompres- 
sion in an overexpanded supersonic nozzle. 

Although many attempts have been made to clarify the aerodynamic mechanism of 
this phenomenon by both theoretical and experimental methods, there does not seem to 
exist any analytical approach pertinent to internal and external flow characteristics from 
which appreciable information is available. The main reason for this deficiency is the 
fact that application of the conventional boundary-layer approximations has been recog- 
nized to be doubtful for this kind of interaction phenomenon. 

Since it seems hopeless to attempt to solve the complete Navier-Stokes equations 
in detail, reference 16 presents an approximate analytical theory of shock-wave — 
boundary -layer interaction without heat transfer derived by modifying the von Karman 
momentum integral method with emphasis being laid on the importance of mass trans- 
port from the outer nearly isentropic stream to the internal dissipative flow through 
mixing. 

In this interaction, it is well-known that the external flow cannot be given a priori 
as a known datum for calculation of the internal dissipative flow. The development of 
dissipative flow itself must be fed back to determine the external flow. This interplay, 
which is now generally recognized as the basis of any valid theoretical treatment, seems 
to have been formulated successfully in reference 16 with the assumption that any state 
of the boundary layer might be predicted by a single parameter. 

For the case of a compressible laminar boundary layer involving separation, the 
validity of this assumption is reasonably verified from the correlation of incompressible 
data through the Illingworth-Stewartson transformation. (See ref. 12.) However, if the 
boundary layer is specified as turbulent, nothing of this kind can be done rigorously 
because of the almost complete lack of theoretical knowledge on the effect of compress- 
ibility, and also because even the empirical information in the incompressible case is 
meager. Despite much discussion on the possibility of considering physical properties 
associated with an incompressible turbulent boundary layer as a one-parameter family, 
it is not yet clear whether this assumption in general is even approximately correct. 
However, if the physical properties do not change too fast, the assumption may be 
adequate. 

Moreover, if heat transfer is involved in the problem, the real physical feature of 
the compressible turbulent boundary layer with a pressure gradient will become much 
more complicated because the energy equation should couple with the momentum equa- 
tion. Since there is neither theoretical nor empirical information available for pre- 
dicting the enthalpy profile across turbulent boundary layers, except for the flat-plate 
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case presented in reference 17, it seems impossible to solve the present problem with- 
out making some assumption concerning the enthalpy profile. 

Introduction of the assumption of a one-parameter family of turbulent boundary 
layers provides a mean temperature within the layer from mass and momentum rela- 
tions, as pointed out in reference 16. However, there is not any rigorous method avail- 
able to examine whether this mean temperature is even approximately consistent with 
the energy equation. 

The purpose of this section is to present an approximate approach to the shock- 
wave — turbulent-boundary-layer interaction with heat transfer so as to emphasize the 
importance of the transport of mass, momentum, and energy in determining the flow pat- 
tern. The present approach is essentially based on the one presented in reference 16 for 
shock-wave — boundary-layer interaction without heat transfer. The basic assumptions 
made in the present approach may be summarized as follows: 

(1) Any state of the turbulent boundary layer with a pressure gradient may be 
predicted as a single-parameter family 

(2) The static pressure across the internal viscous layer is constant and equal to 
the local pressure in the adjacent external flow 

(3) The wall temperature is constant 

(4) The turbulent Prandtl number is constant and equal to unity 

(5) The total enthalpy profile across the boundary layer is assumed to be a function 
of local velocity u only and is given in the form 

H = C P T + | u 2 = (H e - H w )^ + H w (42) 

The final assumption has already been proved to be correct by Van Driest (ref. 17) 
for compressible turbulent boundary layer without pressure gradient if the turbulent 
Prandtl number is unity. In the general case, there is no rigorous proof for equation (42). 
However, it is expected to be approximately valid if the temperature profile within the 
boundary layer does not change too fast. 


Analysis 

To simplify the analysis, a shock wave is assumed to be replaceable by a family of 
isentropic compression waves. This assumption seems to be reasonable because the 
shock wave tends to be diffused and becomes a family of compression waves near the 
boundary layer. A simple sketch of the characteristic flow pattern in the recompression 
region is shown in figure 4. The flow field is assumed to consist of two regions, an 
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External flow direction 
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Figure b.- Flow field in rec empress ion region. 


internal viscous layer of local thickness 6 and an external nearly isentropic flow in 
which the local streamline makes the angle 8 with the x-axis at y = 6. 

The basic equations governing the viscous layer are: 


ax ay 


(43) 


d\X 8U 
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pu — + pv — 
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(44) 

(45) 


Although equation (45) is an energy equation, the solution of this equation has already been 
assumed to be given by equation (42). Therefore, it is not necessary to couple equa- 
tion (45) with equation (44) in the present approach. Thus, equation (45) should be con- 
sidered as an expression for obtaining the heat -transfer rate to the wall. 

Before proceeding to a detailed analytical development, it may be helpful to define 
several physical quantities which will be used in the subsequent argument. Fluxes of 
mass, momentum, and energy across the viscous layer are expressed, respectively, as 


m = 1 

1 pu dy = P e u e (6 - 6*) 

(46) 

‘-j; 

^ pu 2 dy = p e u e 2 (6 - 6* - 0*) 

(47) 

3 -£ 

6 puH dy = He u ~ Hw I + H w m 
) u e 

(48) 
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where equation (42) has been employed in equation (48). From the mass and momentum 
fluxes, mean velocity u within the viscous layer may be defined as 

« _ I _ u e (6 - 6* - e : 

• — 1 i . 

m 6-6 

which provides an important parameter 


(49) 


u _ 5 - 5* - 6 
K u e 6-6* 


(50) 


Crocco and Lees (ref. 16) proposed this parameter as the characteristic one with which 
any state of the viscous layer can be predicted as a single -parameter family. In the 
present approach, the same physical concept as that of Crocco and Lees is to be given to 
this parameter on the basis of the assumption mentioned in the preceding section. A 
mean density p and a mean temperature T within the boundary layer can be defined 
by use of 6 and u together with the equation of state; that is, 


(V _ 

m = pu6 = 


pu5 

RT 


Thus, integration of momentum and energy equations across the boundary layer 
leads to 


cff dm 
dx e dx 
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(51) 
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By use of equation (48), local heat -transfer rate from equation (52) may be expressed as 



- H 


WJ 


/- \ dm ~ dx 

(1 - k ) — m — 

dx dx 


(53) 


On the other hand, geometrical definition of 6 and 6 gives a mass flux relation (see 
fig. 4) 


dm 

dx 



(54) 


which may be reexpressed by introducing the concept of a mixing coefficient, originally 
proposed by Crocco and Lees, as 

17 - ( 55 ) 


31 


L 


Therefore, 


m = k+6 

dx 


( 56 ) 


Since the shock wave is assumed to be replaceable by a family of isentropic waves, 
the pressure in the external flow is given by the Prandtl-Meyer relation: 

dp 


p $ 


M* 2 -! 


(57) 


Moreover, it is convenient to introduce several reduced quantities which are defined as 
u ^ 


w = 


m = ma r 
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^e = 


yw e 


(58) 


By use of these values, equation (57) can be reexpressed as 
dp _ dw e 

P "< P e 

To evaluate several physical properties of compressible turbulent boundary layers 
by use of the incompressible data, the Illingworth-Stewartson transformation (ref. 12) is 
introduced as 

p 0 a o dr ? = P a e d y 

Po a o d £ = Pe a e 
= tv 


o V ~ ~t x 
P o x 


U i = tfj 


Vi = 


and provide a useful correlation between compressible and incompressible flows. The 
results of the transformation may be summarized as follows: 
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P 0 U i d7 7 = P u dyl 
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Po T o 1 




( 59 ) 


where the subscript o denotes stagnation conditions of the compressible external flow, 
to which the corresponding incompressible external flow should be referred. 

However, it must be noted that there is no rigorous proof for the validity of the 
Illingworth-Stewartson transformation in the case of a compressible turbulent boundary 
layer. Because of the almost complete lack of knowledge of compressibility effects on 
turbulent boundary layers with a pressure gradient, and also because no other method is 
available for correlating the compressible turbulent boundary -layer flow with incompress- 
ible flow, the Illingworth-Stewartson transformation is used. 


Several important relations which result from this transformation are summarized 
as follows: 

P p a e , 
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m = p e u e (5 - 6*) = p 0 u ie (6 i - 8**) = m i 


where the subscript i denotes corresponding incompressible conditions. The thickness 
of the boundary layer 6 is obtained by 
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If local stagnation temperature in the viscous layer is denoted by T s , then 
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and therefore, 


T _ T T s _ H y - 1 w 2 
T 0 T s T 0 H e 2 

which can be reexpressed by use of equation (42) as 


= X + (1 - X) 


JL . 1 - 1 w 2/JL 


u c 




With the aid of equation (59), the local thickness of the compressible boundary layer is 
written as 

6 = ^ H 6 * ■ (I ■ x)6 ‘* ■ ^ We2 ^ ■ h * ■ Si *l 


where 
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With this value, the mean temperature may be obtained in a form 
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Thus, the mass flux relation (eqs. (58)) is reduced to 

- (1 - X) F 2 + 1 - w e 2j = ypw e 6 
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Differentiation of equation (60) with respect to x and elimination of d6/dx by use of 
equation (56) leads to 
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where 


= Fj - (1 - X) F 2 + 1 - 2^-=-^ w e 2 


(62) 


On the other hand, with the aid of equation (49), the momentum equation (eq. (51)) 
may be reduced to 
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Elimination of dw e /d(i Q g m) from equations (61) and (63) yields for (p ¥= 1 
d k B 


dw e A 


where 
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We 2 )«l-^ 


^ w e 2 )| + « (1 - (‘ - H 1 w e 2 ) 


[f x - (1 - X) F^j (1 - 0 + K<t>) - k(1 - k)( 1 - 0) 


(63) 


(64) 


(65) 


Hi . (1 . JO 5!£Tl 

_dK ( } d/Jj 


B = k( 1 - /c) (1 - 0)^1 - We 2 ^ 21 ^- w e 2 - l) + k [F x - (1 - *)F 2 ] 

x |Vw e 2 (l - k)( 1 - 0) + 0(1 - k)^ 1 - Wg 2 ) 

+ (» - ^ w e 2 )| - «[ F 1 - d - « p 2]j 

Equation (65) is a fundamental equation associated with the rate of change in the state of 
the internal viscous layer with external conditions. As is seen in the equation, the essen- 
tial characteristic of the interaction phenomenon is nonlinear. If 0 is equal to unity, 
the equation can be reduced to a simple form, 
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(66) 


d k _ ~ (1 ' X)F 2 

dw e w e( 1 ~ Z Y^ w e 2 ) 

which is independent of mixing coefficient. 

Crocco and Lees (ref. 16) gave a detailed discussion on the solution of equation (65) 
in the case where X = 1.0 and 0 = 0 and pointed out the existence of a singular point 
which is a saddle point of the solution from a mathematical viewpoint. Physically 
speaking, this singular point corresponds to the sonic point of a convergent-divergent 
nozzle. However, if 0 is equal to unity, the singular point vanishes, as is seen in 
equation (66). 

Evaluation of and F 2 is of interest. Crocco and Lees used an empirical 
relation for Fj, obtained from the data given by reference 18. However, in the present 
approach, F^ and F 2 are estimated from approximate analytical results obtained for 
the fully developed turbulent boundary layer on a flat plate. If a velocity profile for an 
incompressible turbulent boundary layer with zero pressure gradient is assumed to be 
given by a power-law distribution such as 

u j = m 1/n 
u ie \^i J 

then, the displacement thickness and the momentum thickness are obtained, respec- 
tively, as 



6 i 1 + n 


^1 _ n 

«i “ (1 + n)(2 + n) 


where n indicates a positive constant. For the flat plate, it is usually assumed that 
n = 7. With these expressions, k, Fj, and F 2 are given in forms: 


n + 1 
+ 2 




F 





(67) 


F 2 = 


2 + n 
n(l + n)^ 


Therefore, elimination of n from these equations leads to 
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( 68 ) 


Fi = g(l ^ 

1 2k - 1 


F 2 = 


k(2k - 1) 


Although equations (68) were derived from a relation available for a flat plate, it is 
capable of representing a universal relation associated with the general state of the 
incompressible turbulent boundary layer, as long as the assumption of a one-parameter 
family of boundary-layer state is valid. The reason for this statement is as follows: 

From the assumption of a one -parameter family for the boundary-layer state, it is 
recognized that there should exist a unique relation between Fj and k which can 
predict any state of the boundary layer. In other words, Fj should be a function of k 
only. Since flat-plate flow having k = /c p and Fj = F^p for a specified value of n 
must be one of the general states, these values have to coincide with those of the unique 
solution. This condition, in turn, means that the parameter n in equations (67) may be 
regarded as a parametric variable indicating various states of the boundary- layer flow. 
Therefore, equations (68) which are obtained by eliminating the parametric variable can be 
regarded as a representative form of the unique relation between Fj and k. The same 
reasoning is true for F£. 


Evaluation of </> in equation (64) is another point of interest, 
it is easily shown that 


From equations (59) 


C f = 


p e p o a e^/ u ie 


p o p e a 
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On the other hand, from the definition of the mixing coefficient 


dm 

dx 


= Kp u e = ££ = KjP u ie 

J -- e e dx 1 o ie 


P e a e 

P Q a O 


Therefore, 


K = 


p ePo K 

p p n 

F o^e 


(69) 


and 


0 = 


Cf 


C f,i 


2K(1 - k) 2Kj(l - Ki) 


= <Pi 


(70) 


Equation (70) clearly indicates that <fi does not change through the Illingworth- 
Stewartson transformation. For the flat-plate flow, the incompressible mass flux is 
given by 


L 
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ffi i = P 0 u ie( 5 i " 6 i1 = Po u ie YT^ 6 i 

Differentiation of rrij with respect to £ leads to 


Ki = 


n d6 i 


1 1 + n d£ 

On the other hand, the skin-friction coefficient is given by 
dft* 


C t . = 2 


2n 


d5 4 


d£ (1 + n)(2 + n) d£ 


( 71 ) 


Hence, 


and 


2lA - 1 = i _ K . 

2Ki ' 2 + n Kl 




for a flat plate. Unfortunately, for general states of boundary-layer flow other than the 
flat plate, there does not exist any rigorous theory for evaluating 4> Moreover, even 
the empirical information is meager. For this reason, it is assumed in the present anal- 
ysis that (p is equal to unity everywhere in the recompression region. 

The value of the mixing coefficient remains to be discussed. The mixing coeffi- 
cient of a compressible family is correlated with that of the incompressible one through 
equation (69). It is easily deduced from equation (71) that the local mixing coefficient 
for attached boundary-layer flow is, in general, a function of the local Reynolds number. 
However, even if such a formulation had been made, it is still doubtful whether it could 
give a correct answer that would include the effects of compressibility. When it is con- 
sidered that the mass-flow variation has only a little effect on the turbulent values and 
that insufficient empirical data are available on the spread of turbulence, it seems justi- 
fiable to introduce the approximation that the mixing coefficient is constant throughout 
the recompression region. 

From the data presented in reference 19, Crocco and Lees estimated that 
K = 0.03 near the separation point of the attached turbulent boundary layer. A rough 
evaluation of the mixing coefficient by use of equation (29) indicates that K « in 

the jet mixing region. Of course, this value may vary depending upon the definition of 
the upper jet boundary. Nevertheless, it is more than four times the Crocco and Lees' 
value for X = 1.0. On the other hand, the flat-plate value for a Reynolds number of 10^ 
indicates that K « 6 X 10"®, which is only one-fifth of Crocco and Lees' value. When 
the fact that the actual value of the mixing coefficient in the recompression region should 
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vary from K » at the entrance to K « 6 x 10 at the exit of this region is con- 

sidered, the value proposed by Crocco and Lees is realistic as a mean value. Therefore, 
the mixing coefficient is, for the present, taken as K = 0.03, where this value should be 
regarded as a realistic mean value in the recompression region. However, it should be 
noted that the value of the mixing coefficient has a strong influence on the thickness of 
the boundary layer and the associated local heat -transfer rate in the recompression 
region. 


Since the viscous layer downstream of the recompression region has been assumed 
to be a fully developed turbulent boundary layer on a flat plate, integration of equation (65) 
is made upstream starting from the exit point of the recompression region with the known 
value of w e and k appropriate to the flat-plate flow conditions. From the assumption 
of isentropic recompression, the Mach number of the external flow downstream of the 
recompression region must equal the oncoming value. Therefore the starting values 
are given by 


(w P ) = w Q „ 
v e; Start e >P 


^ Start K p 9 






V 2 




J 


(72) 


For the case of <p = 1 , equation (66) can be integrated analytically in a form 


— — 1 / 

(1 - k) 1+X (2k - 1 + X) 1+X = C I 


JLl w 2 
2 e 


Wp 


(73) 


where C is an integral constant to be determined from equations (72). 


Other Relations Associated With Recompression Region 

Since the solution of the fundamental equation (eq. (65) or eq. (66)) has been 
obtained, other quantities concerning the recompression region can be estimated. For 
this purpose it will be convenient to introduce nondimensional quantities defined as 
m -\ 


m = 


m r 


s 5 

6 = ^ > 


X = -■ 


( 74 ) 
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where the subscript p denotes conditions at the exit of the recompression region and 
the origin of the x-axis has been taken at the exit point. 

The equation for mass flux is given from equation (63) for <p * 1 as 


dm _ m 

dw e w e (l - /r)( 1 - 0) 


e 


d/c 

dw e 


/c[f x - (1 - X)Fj| 

yw e 0 e 


In the case 0 = 1, this equation may be reduced to 


dm ^2 
dw e Ei 


(75) 


(76) 


where 



and where P has been given by equation (62). The boundary condition is given by 

m = l ( w e = w e,p) ( 77 ) 

By use of the mass flux m r thus obtained at the entrance of recompression region, 
the value of which was unknown in the earlier development, can be determined through 
continuity of mass flux between the jet mixing region and the recompression region. 

From equation (60), there results 

/5_\ _ K r p r 

VmJr yPjWgj (78) 

where the subscript r denotes conditions at the entrance of the recompression region 
(or exit of the jet mixing region). On the other hand, the mass flux at the exit of the jet 
mixing region has been given by equation (28), and the thickness of the mixing layer at 
the same station is given by the equation 


6 r = 



L 

a l\|%e^L 




Therefore, 

(±\ / 6 \ J(gu)-J(?i) 

U/ r {mz 0 ) r Pj Vo( C U ' ? i) 


(79) 
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Thus, from equations (78) and (79), the following relation which determines 
uniquely is obtained: 


J (^u)~ J ( £j) _ K r[fl >r - (1 - X)F 2>r + 1 - w e>1 


?u-?i 


1 . Y— 1 


2 W e,l" 

Thickness of the boundary layer is estimated from equation (60); that is, 

6 Pp w e 


6 - -2- - ISiE JL Z.S 

6 p P w e K p Pp 


( 80 ) 


which is further rewritten by use of isentropic relations in a form 


- w p n 

6 = e ?P 
w c 



(81) 


The equation for the width of the re compression region can be obtained from the defini- 
tion of the mixing coefficient. Rewriting equation (55) yields 


dm _ 5 p 
dx m p 


P e u e a o K 


which can be reexpressed, with the aid of equation (60), as 


dx 

dw e 


w e,p 1 ~ 1 T ~ V 


^ K P P P 


w e 


11 - 1 ~ 1 w 2 
' 2 e,p 
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1 -Z_^i 


w c 


dm 

dw„ 


(82) 


where the boundary condition is given as 

x = 0 ( w e = w e>p ) 

The local heat -transfer rate can be obtained from equation (53) as 


l w 


- y z 1 


w, 


Poo U co(H e “ H w ) 


e,p 


k p p p 


dw e 


dm 

dw e 

dwg 

dx 


(83) 


(84) 


The mean heat -transfer rate q R in the recompression region is defined by the equation 

>0 


q R 


" x r J Xt qw 


dx 
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where the negative sign of x r is due to the coordinates system used for the recompres- 
sion region. (See fig. 4.) By use of nondimensional expressions, the preceding equation 
can be reduced to 


q R 


P<=o U °°(**e ~ H w ) x r 


« dW£ 

x r ^w e>p q ’ dw e 


(85) 


where 

C = — 

q,R ” Pco U =°( H e - H w) 

which has already been given by equation (84). The mean heat -transfer coefficient in the 
recompression region C q R is defined by the equation 

C = 

q »^ P^UooHco 

HEAT -TRANSFER REDUCTION IN THE REGION 
DOWNSTREAM OF RE COMPRESSION 


Since the coolant flow is assumed to be injected into the dead-air region at a low 
flow rate through a rearward-facing slot, complete mixing occurs between the coolant 
and the oncoming viscous flow before the exit of the recompression region is approached 
so that the viscous layer downstream of the recompression region is an ordinary turbulent 
boundary layer on a flat plate (fig. 5). Therefore, the predominant factor in reducing the 
heat-transfer rate in this region will be the growth of boundary -layer thickness. 

Since the recompression has been assumed to be made through a family of isen- 
tropic compression waves, the Mach number of the external flow must be the oncoming 
value. In contrast to the jet mixing region upstream, this characteristic makes the sub- 
sequent analysis fairly simple because the external flow Mach number remains unchanged 
with coolant flow rate. 
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Figure 5.- Basic turbulent boundary layer along 
an infinite flat plate. 



For fully developed compressible turbulent boundary-layer flows without pressure 
gradient (flat plate), Van Driest (ref. 17) presented an analytical approach whose results 
are shown to be in good agreement with existing experimental data. Therefore, subse- 
quent development will use the analytical results presented by Van Driest. It must be 
noted that although coolant flow rate does not appear explicitly in the subsequent analysis, 
its effect is implicitly involved in the development because the behavior of the attached 
boundary layer is seriously influenced by the upstream conditions. 

The local skin-friction coefficient Cf ^ and the energy transfer coefficient Cjj )OQ 
obtained by Van Driest for a fully developed turbulent boundary layer on a flat plate may 
be summarized as follows: 


(sin -1 Q! + sin -1 ^) = log 


N l C f,°° 1/2 ' Tw ' 
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and where x denotes a distance measured from leading edge of the flat plate where the 
thickness of the boundary layer is zero. It should be noted that these relations are 
obtained under the assumption of a turbulent Prandtl number equal to unity. Thus, the 
local thickness of the boundary layer is given by 


6 = 


_ K(l + N 2 - Nl^ /T,y /2 C 1/2 

1/2 It 


C f "x 

1,00 


where 


K = 


(91) 


(92) 


4.15\/2 log jq e 

Therefore, equation (91) can be simplified, by use of equation (90), to a form 
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= -!_ xC , 1/2 


3.6 


L ,0O 


(93) 


On the other hand, from equations (86) and (89), a relation between the local skin- 
friction coefficient and x can be obtained as 
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exp 
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Substitution of equation (93) into equation (94) and elimination of x then leads to 
l+2o> 

I z_ . 1 /O I 
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(94) 
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UqqD 


(96) 


Thus, if the Mach number of the external flow Moo and the temperature ratio T^T^ 
are given, solutions for the basic turbulent boundary -layer flow (Van Driest’ s solution) 
become known in detail. 

For the present problem, a fully developed turbulent boundary layer has been 
assumed to develop downstream starting from the exit of the recompression region with 
a finite thickness 6 p . Based on the discussion in the previous paragraphs, it seems 
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reasonable to assume that characteristics of the actual boundary layer are the same as 
those of the basic boundary -layer flow, in the region downstream of a point where the 
thickness of the basic boundary layer is equal to 5p. Furthermore, 6p/L is known 
from the solutions of the jet mixing region and recompression region to be 

fp _ _l_£r 
L 5 r L 


Therefore, the corresponding station Xp and the skin-friction coefficient Cf w p in 
the basic boundary layer are found from equations (93) and (95) to be 


l+2o> 



If the coolant flow rate increases, the other oncoming flow conditions remaining 
unchanged, 6p will increase because of the increase of mass flux across the boundary 
layer. This increase, in turn, tends to shift the starting point of the actual boundary 
layer downstream when it is transferred to the basic boundary-layer contour. (See 
fig. 5.) This statement follows, since, as is seen in equations (97) and (98), the increase 
in 6 r results in a decrease in Cf >00> p and, consequently, an increase in Xp*. More- 
over, since the skin-friction coefficient decreases, the surface heat input will to a cer- 
tain extent decrease. 


In order to examine this effect more quantitatively, the local effectiveness of cool- 
ant injection r/f associated with the fully developed turbulent boundary -layer region is 
defined as 


aOO OO 
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( 100 ) 
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Since the local heat -transfer rate q w (x) associated with the basic boundary-layer flow 
is kept unchanged irrespective of the coolant flow rate, the numerator on the right-hand 
side of equation (100) may be evaluated as 


C°° r°° P x p _ 

1 Vo^-l q w ' J ’.W* 

A p,o A p A p,o 


where q,..(x) has been given by equation (87) and the subscript o denotes the no- 
w 

coolant-flow condition. Thus, by introducing nondimensional expressions and replacing 
x by Cf >oo as an integral variable, equation (100) may be reduced to 


where 



( 101 ) 


x exp 


f~~) (sin -:1 Q! + sin“l/3) - 0.944 


( 102 ) 


The overall effectiveness of coolant injection may be defined by an equation 

f(Vo - %)«* 


Tjrri — ■ 


Pc u c H ooD 


(103) 


where the origin of the x-axis is taken at the slot exit and q w denotes local heat- 
transfer rate to the wall. Equation (103) can be reexpressed as 
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where L and l denote the length of the jet-mixing and recompression regions, respec- 
tively, and the subscript o indicates the no -coolant -flow condition, and where subscripts 
j, R, and f denote conditions in regions of jet mixing, recompression, and fully devel- 
oped turbulent boundary layer, respectively. 
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RESULTS AND DISCUSSION 


Jet Mixing Region 

Numerical calculations for determination of the flow field in the jet mixing region 
can be carried out by the following procedure. First, the Mach number of the oncoming 
flow Moo, the ratio of wall enthalpy to stagnation enthalpy of oncoming flow H w /Hoo, and 
the ratio of stagnation enthalpy of coolant flow to wall enthalpy H c /H w are specified, 
and a flow deflection angle 9 of the external isentropic flow through a Prandtl-Meyer 
fan is assumed. These values permit calculation of the Mach number of the external 
flow Mi in the jet mixing region, and also determination of the shock- wave strength in 
the recompression region. It must be noted that the assumption of 9 specifies the 
pressure p^ in the jet mixing region. Since the strength of the shock wave is known, 
any solution associated with the recompression region can be considered to have already 
been determined at this stage of calculation. 


Since Mi is known, velocity and enthalpy profiles are calculated from equa- 
tions (16) and (17), respectively. These profiles, in turn, permit determination of 
by use of equation (80). If u c * is assumed to be zero (no coolant flow) at this flow 
deflection angle, then the parameter Z becomes known from equations (24) and (25), 
because Thus, by changing the value of 6 (or Mi) but keeping u c * = 0, a 

relation between p-. and Z /= — ^ -~ \ can be obtained. This relation indicates the 


usual base pressure associated with a rearward-facing step, which should be the mini- 
mum value of pressure in the jet mixing region for the case of film cooling. Since 
oncoming flow conditions and slot geometry have been given or are known a priori, the 
value of Z is known and, consequently, 9 or Pi is known from the relation just 
obtained. The data (M 0 , p 0 , Z) thus obtained for the no -coolant -flow case are con- 
sidered to be starting conditions for the coolant-flow case. 


In the case of coolant mass injection, essentially the same procedure as used for 
the no-coolant-flow case is applied to obtain except that Z is kept constant irre- 
spective of coolant flow rate, since the oncoming flow is not influenced by the coolant 
mass injection. For this purpose, a new flow deflection angle 9, smaller than 0 O , is 
specified and the coolant flow velocity is assumed. However, at this stage of calculation, 
it is not yet known whether the assumed value of u c * is compatible with the specified 
value of 9 or Mj. 

Since has been specified, it is easily shown from isentropic relations that 

the static enthalpy of coolant at the slot exit is given by the equation 
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(104) 


Therefore, an assumed local coolant flow rate can be obtained from the assumed value 
of u c * through the equation 



Moreover, 6 ^ * JT) and 0^*/D are obtained from equations (26) so that and J(? m ) 
are known from equations (23) and (25). 

Once and J(£ m ) have been determined in this way, the length of the mixing 
region L/D is known from equation (21) and, consequently, p c * u c * can be evaluated 
again from equation (19). However, this resultant coolant flow rate does not always coin- 
cide with the value initially assumed because, physically speaking, the solution to u c * 
must be unique for a specified value of 0. Therefore, it is necessary to repeat the cal- 
culation, the assumed value of u c * being adjusted, until the resultant coolant flow rate 
coincides with the assumed one. The coolant flow velocity and coolant flow rate obtained 
in this way are unique for the specified values of 0 and Z. 

Figure 6 shows a typical example of velocity and enthalpy profiles across the mixing 
layer for M*, = 3.0, 0 = 9° (Mj = 3.514), and = 0.3571. Figure 7 presents the 

variation of pressure in the jet mixing region with Z for Moo = 2.0, H w /Hoo = 1.0 and 
u c * = 0, where the oncoming boundary layer is assumed to be laminar (n = 2). In the fig- 
ure are also plotted the experimental results obtained by Chapman, Wimbrow, and Kester 
(ref. 15) for M^ = 2. However, since the value of <£j expressed by equation (31) has 
been determined to obtain good agreement of base pressure between the present approach 
and the data of reference 15 at M*, = 2, the good agreement indicated in figure 7 should 
have been a natural result. Because of this circumstance, base pressure has been calcu- 
lated for Moo = 3 to confirm the present approach, and the results are shown in figure 8, 
together with the data of reference 15 for comparison. As is seen in figure 8, agreement 
between theory and experiment is fairly good; thus, the present approach is confirmed. 

Figure 9 shows several examples of pressure variation in the jet mixing region with 
coolant flow rate, where Z is used as a parameter. The pressure increases with the 
increase of the coolant flow rate. This characteristic may be understood by the following 
consideration. If parametric conditions are fixed, it is easily shown that the minimum 
pressure attainable in the dead-air region is the base pressure corresponding to u c * = 0, 
below which fluid in the dead-air region will be sucked into the coolant slot. Therefore, 
in order that the coolant mass be injected, the pressure of coolant flow at the slot exit 
must be higher than the base pressure without coolant flow. This condition, in turn, 
increases the pressure in the dead-air region, since coolant flow is subsonic at the slot 
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Present theory 







Figure 9-- Variation of pressure in jet mixing region 
with coolant flow rate. = 3-0; H w/Hw = 0 . 3571 ; 

H c /H w ” 1*°; Np r ,t = 1-0* 


exit. The increase of dead-air pressure decreases the strength of the Prandtl-Meyer 
fan. More intuitively, it might be deduced that the separated boundary layer, which acts 
as though it were an impermeable membrane, is so pushed upward, through a higher 
pressure inside which is supported by the coolant flow, as to reach a new equilibrium 
position and, therefore, the flow deflection angle of the external stream decreases to 
increase the pressure downstream of the Prandtl-Meyer fan. 


Figure 10 presents a variation of the length of the separated region (length of 
mixing region) with coolant flow rate for Moo =3.0, H w /H OO — 1.0, and Z = 2.8 X 10“1 
As is seen in the figure, the length increases with an increase of coolant flow and this 
characteristic supports the intuitive deduction of a pressure increase, due to the coolant 
flow, mentioned in the last paragraph. 


In order to confirm further the validity of the present approach, an experimental 
curve for effective length of film cooling obtained by Goldstein, Eckert, Tsou, and 
Haji-Sheikh (ref. 4) is shown in figure 10 for comparison. However, it must be remarked 
that there is a difference in the physical meaning of L/D between theory and experi- 
ment indicated in figure 10. In the experiment L is defined as the length measured 
from the slot exit for which the effectiveness of film cooling 77^ defined by the equa- 
tion 


^eff 


• w 


t* 

j. « 


T w,2 " T r,2 


is assumed to be kept as unity, where T w is the adiabatic wall temperature and T w 2 j 
its value at the point of injection. Here, T' is the wall temperature for isoenergetic 
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Figure 10.- Length of separated region. =3*0; 

Hv/H^ = 1.0; H c /H w = 1.0; Z = 0.0282; 

X = 1.0; N pr ^ t = 1.0. 

injection (H c = Hoo) at the same blowing rate, and T^, ^ indicates the value at the point 
of injection. Thus, the experimental length L/D in figure 10 has the physical meaning 
that the wall temperature T w changes abruptly to approach T^, downstream of the 
point where x = L. However, in spite of the physical meaning of L/D being different 
between the present theory and the experiment, the comparison shown in figure 10 seems 
to be reasonable because the change of wall temperature clearly indicates an abrupt 
change of surface heating rate which will take place at the end of the jet mixing region, 
because the mixing between oncoming flow and coolant flow has almost finished before 
the recompression region is approached. This reasoning can be confirmed by observa- 
tion of schlieren photographs of the shock-wave location presented in reference 4. Thus, 
the experimental length L/D might be interpreted as the length of separated region. 

Goldstein, Eckert, Tsou, and Haji-Sheikh plotted the value of ri e ^ against a 
parameter \ defined by the equation 



and obtained \ = 23.5 for ri e ^ =1 by means of extrapolation of the mean data curve 
in the case of small coolant injection (M % 0.12). However, detailed examination of the 
experimental data (fig. 4 of ref. 4) clearly reveals that f = 23.5 is too high; the correct 
value should be near \ = 20. Since a choice of | = 20 seems to be more realistic from 
the experimental viewpoint, an empirical curve defined by the equation 

L = 20M°" 4 (106) 
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is also shown in figure 10 for comparison. As is seen in the figure, agreement between 
present theory and the empirical data given by equation (106) seems to be fairly good. 

Figures 11 and 12 show the effects of Z and wall enthalpy on the length of the 
separated region. In general, the length of the separated region primarily depends upon 
the flow deflection angle which the external stream makes with the wall, and this condi- 
tion leads to the intuitive deduction that the weaker the Prandtl-Meyer expansion at the 



Figure 11.- Length of separated region. H*, = 3-0; 
Hw/H^ = 0.3571; Hc/H^ = 1.0} X = 1.0; 

N Pr,t = 



Figure 12.- Length of separated region. =3-0; 

H c /H w = 1.0; Z = 0.0^2; X = 1.0; N pr ^ t = 1.0. 
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edge of the plate lip, the longer the separated region becomes. Since the strength of the 
Prandtl-Meyer fan decreases with the increase of either Z or M as is seen in fig- 
ure 9, the results shown in figure 11 seem to be quite reasonable. 

It is generally recognized that the thickness of an attached boundary layer becomes, 
smaller for a lower wall enthalpy. Detailed examination of velocity and enthalpy profiles 
across the jet mixing region reveals that the same is true for turbulent jet mixing. In 
fact, the value of the second term on the right-hand side of equation (21) decreases as the 
wall enthalpy decreases and, consequently, the length of the separated region increases. 
However, it is easily shown from a comparison of figure 12 with figure 11 that the length 
of the separated region is less sensitive to wall enthalpy than to the parameter Z. The 
same trend has already been pointed out experimentally in reference 4 for the case of 
coolant injection through a rearward-facing slot; thus, further experimental evidence is 
given of the validity of the present approach. 

Figure 13 presents the variation of the mean heat -transfer rate with Z for no 
coolant flow. It is well-known that the heat -transfer rate for an attached boundary layer 
depends primarily upon the rate of energy transport into the boundary layer from the 
external stream, the difference between adiabatic wall temperature and actual wall tem- 
perature, and the boundary -layer thickness. Therefore, if the wall temperature and the 
stagnation enthalpy of the external isentropic stream are assumed to be fixed, then the 
local heat -transfer rate depends mainly upon the local Mach number of the external 
stream. This condition follows because the local thickness of the boundary layer and the 
mixing coefficient in general do not change so abruptly even under the influence of pres- 
sure gradient, except for some particular cases such as flow separation. Inasmuch as 
the local rate of energy transport is proportional to the local mass flow rate of the exter- 
nal flow, which decreases with an increase of local Mach number if the flow is supersonic, 
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Figure 13.- Mean heat-transfer rate in jet mixing 
region. M* = 3-0; E w /'B m = 0.3571; M = 0; 
%r,t = 1 -°- 
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the local heat-transfer rate is expected to decrease considerably as the local Mach num- 
ber increases. This fact will be recognized more readily by considering the physical 
evidence that the maximum heating rate is encountered in the region of a sonic throat of 
a convergent-divergent supersonic nozzle. 

In the case of an attached turbulent boundary layer, the local rate of energy trans- 
port increases considerably when compared with that of an attached laminar boundary 
layer, because the local mixing coefficient becomes of the order of ten times the laminar 
value for the same Reynolds numbers. Moreover, the growth of boundary-layer thick- 
ness has less effect on the local heating rate than the case with laminar flow, because 
q w oc 6 “1/4 for turbulent flow whereas q w «: 6“* for laminar flow. These unfavorable 
characteristics seem, therefore, to be largely responsible for the marked difference in 
heating rate between turbulent and laminar boundary-layer flows. 

In the case of turbulent jet mixing which follows flow separation at a rearward- 
facing step, there seems to exist at least two favorable factors, where heat transfer is 
concerned. The one is that the external flow Mach number increases to a certain extent 
from its oncoming value through the Prandtl-Meyer expansion. Another is that the effec- 
tive thickness of the viscous layer increases abruptly to the extent of the step height, 
although this effect may not be so efficient for turbulent boundary layers. In spite of 
these favorable factors, however, the turbulent jet mixing region has a marked unfavorable 
factor in that the local mixing coefficient increases considerably when compared with its 
value for attached boundary layer. This unfavorable factor seems to overcome the favor- 
able factors and results in a net increase of heating rate in this region. 

Figure 13 indicates that the mean heat -transfer rate in the jet mixing region 
increases with an increase in Z. This increase can be interpreted to mean that the 
heating rate increases with either an increase of oncoming boundary-layer thickness or a 
decrease of coolant slot height. The effect of coolant slot height can be understood rather 
easily. However, the effect of oncoming boundary-layer thickness given by the present 
approach may seem to be inconsistent with the general characteristics of the heat-transfer 
rate. 

As has already been mentioned, strength of the Prandtl-Meyer expansion fan depends 
upon the value of Z if the oncoming flow Mach number is fixed and the coolant flow does 
not exist. Since the external flow Mach number M^ decreases with an increase in Z, 
the rate of energy transport and the static temperature at the edge of the mixing layer 
also increase. The favorable factor of increasing thickness of the oncoming boundary 
layer is not very effective and the unfavorable factors predominate and result in a net 
increase of heating rate in the jet mixing region. Thus, the qualitative trend of the result 
shown in figure 13 is reasonable. 
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However, it is necessary to determine whether the quantitative magnitude of the 
mean heating rate is reasonable. Unfortunately, there exist neither theoretical results 
nor empirical data to compare with the present result. For this reason, the mean heat- 
transfer rate is approximated by use of a simple analogy, based on the result for an 
attached turbulent boundary layer. Detailed examination of turbulent boundary-layer 
characteristics on a flat plate reveals that the local heat -transfer rate is primarily pro- 
portional to the local mixing coefficient and the thickness of boundary layer in a form 



Although the aerodynamic mechanism of separated flow may be essentially different from 
that of the attached flow, it is assumed for convenience that the heat -transfer rate asso- 
ciated with the separated flow can be expressed by this relation. Since the external-flow 
Mach number in the jet mixing region is different from its oncoming value, the effect of 
mass flow rate must be taken into consideration. Thus, the ratio of heat -transfer rate 
between the separated flow and the associated attached oncoming flow at the point of sepa- 
ration is given by the equation: 

Vs Ks/M 1/4 Pl u l 

Va Ka W P ” U ” 

where the subscripts s and a denote conditions appropriate to separated and attached 
flows, respectively, and the subscripts 1 and °° indicate conditions of external flow 
pertinent to jet mixing and oncoming regions. For simplicity, 6 a is assumed to be 
given by the thickness of the oncoming boundary layer 5oo, and 5 S is assumed to be 
given by the sum of 5,*, and the coolant slot height. Therefore, 


£a 

D 


5qq ^ 

D 


5Z 


5g = Soo + D 

Under the conditions 
Moo = 3.0 

= 0.3571 

Hoo 

Z = 0.042 
N Re,D = 1q4 
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the oncoming boundary -layer thickness is estimated as 

= 0.21 

D 

The present approach gives as the value of mixing coefficient for separated flow 
0.139 


K s = 


X 


which seems to be nearly constant irrespective of the given conditions. 

On the other hand, the mixing coefficient and the associated heat -transfer rate for 
attached oncoming flow are evaluated at 5oo/D = 0.21 by use of flat-plate solutions given 
by Van Driest (ref. 17), which are obtained, respectively, as 

K a = 0.017 
qw * a -= 1.49X10-3 

p UcoHco 

O O 

Therefore, the estimated heat -transfer rate for the associated separated flow can be 
obtained for X = 1 as 


= 5.54 x 10" 3 

P^UooHoo 

Similarly, the estimated heat -transfer rate can be obtained for other values of Z and 
the results are given in table I together with the present results of mean heating rate 
within the mixing region for comparison. As can be seen in the table, the order of mag- 
nitude of the present results is comparable with the estimated values, except that the 
present results are slightly larger than the estimated ones. However, in view of the fact 
that the estimated value corresponds to the local heating rate just downstream of 


TABLE I.- COMPARISON OF THE PRESENT RESULTS 

WITH ESTIMATED HEAT-TRANSFER RATE 

[M = 3.0, Hw/Hoo = 0.3571, K s = 0.139/X, 

M = 0, N Pr>t = 1.0] 


7 

w,s/£ 

*oo U °°^°o 

^w,s/lw,a 


Present 

Estimated 

0.042 

5.99 x 10 3 

5.54 X 10 3 

3.72 

.082 

7.78 

6.92 

5.40 

.101 

8.38 

7.32 

l 

6.00 


I 


separation point whereas the present result indicates a mean heating rate with the mixing 
region, small quantitative differences between these results are not important because 
the estimated value itself is only an approximate result. In any event, it seems that the 
present results will not differ much from the exact heat -transfer rate in the jet mixing 
region, if X is reasonable. 

It is seen that the heating rate increases considerably in the separated region, when 
compared with that of the attached flow. This fact has already been pointed out by 
Chapman (ref. 20). According to his estimate, the heating rate for separated turbulent 
flow becomes nearly six times the attached flow value at Moo = 0 and this ratio 
decreases to 2.8 times that in the corresponding attached boundary layer at Moo = 1.6; 
thus the effect of Mach number is strong for turbulent flow. Chapman (ref. 21) further 
deduced that although this characteristic may be in sharp contradiction to the analogous 
comparison for the laminar case, the ratio will decrease considerably as the Mach num- 
ber increases sufficiently. Since Chapman's evaluation is pertinent to the case of zero 
thickness of the oncoming boundary layer (Z = 0), a direct comparison between the pres- 
ent result and Chapman's data cannot be made. However, when it is considered that the 
ratio decreases with a decrease of Z, the present value of this ratio for Z = 0.042 
seems to support Chapman's deduction. 

However, it must be noted that there still remains a slight uncertainty in the quanti- 
tative value of the mean heat -transfer rate obtained in the present approach. Since the 
mean heating rate in the jet mixing region depends strongly upon the value of o^\jNj^ e ^ 
as is seen in equation (38), the choice of this value will modify the quantitative value of 
the mean heat -transfer rate. Although it has been evaluated by use of equation (15), there 
does not seem to exist in the present approach any rigorous proof that X = 1.0. In this 
sense, the present approach could be improved with better information on the flow prop- 
erties of compressible turbulent jet mixing. 

Figure 14 shows the effect of wall temperature on heating rate. As can be seen in 
figure 14, the mean heat-transfer rate without coolant flow decreases with an increase of 
wall enthalpy, and this trend is quite reasonable. 

Figure 15 presents the variation of the mean heat -transfer rate in the jet mixing 
region with the coolant flow rate, and the effect of wall enthalpy is examined. Figure 15 
clearly indicates that the coolant mass injection is effective in reducing the heating rate 
in this region, but the effectiveness itself seems to be small compared with the laminar 
case at the same blowing rate. This condition occurs because the energy input entrained 
from the external flow is much larger for turbulent mixing than for the laminar case. It 
can be seen that coolant mass injection becomes increasingly ineffective as the wall 
enthalpy decreases. This condition may result because the mixing layer becomes thin. 
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Figure 16 shows the effect of the parameter Z on the reduction of the mean heat- 

transfer rate. In figure 16, it can be seen that q ,/q . is rather insensitive to Z. 

Since the external flow Mach number varies less with coolant flow rate as Z grows, 

the rate of energy input from the external stream to the internal viscous layer becomes 

small for large values of Z. This characteristic seems to increase the effectiveness 

of coolant mass as a heat sink and, consequently, Qw,j/^w,j,o becomes insensitive to Z 

even if q w i n itself is large for a large value of Z. 

W ?J? U 



Figure l6.- Variation of mean heat-transfer rate in jet mixing 
region with coolant flow rate. Moo = 3-0; H w /Hco = 0-3371; 

H c /H w = 1*0; X = 1 • 0 j Npr^t = 1*0* 

In order to demonstrate the effect of mixing coefficient on flow characteristics in 
the jet mixing region, the length of the separated region and the reduction of the mean 
heat -transfer rate with coolant flow rate were calculated for X = 2 and the results are 
shown, as examples, in figures 17 and 18, respectively; the results for X = 1 are also 
shown for comparison. As is seen in figure 17, the length of the separated region 
increases as the mixing coefficient decreases (large value of X). This condition is 
clearly due to the fact that the small mixing coefficient gives a small mixing-layer thick- 
ness and, consequently, elongates the separated region. However, as has been mentioned 
previously, the length of the separated region depends primarily upon the flow deflection 
angle which the external isentropic stream makes with the wall; as a result, the magni- 
tude of the mixing coefficient has a relatively small influence on the length of the sepa- 
rated region. 

On the other hand, the heat -transfer rate in this region is strongly influenced by the 
magnitude of the mixing coefficient. For instance, figure 13 shows that the mean heating 
rate for no coolant flow is inversely proportional to X and, hence, directly proportional 
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to the mixing coefficient. In the case of coolant flow, the mixing coefficient strongly 
influences the reduction of the heating rate. (See fig. 18.) However, this result is quite 
natural because energy entrained in the mixing layer from the external stream decreases 
in proportion to the decrease of the mixing coefficient, whereas the effect of coolant mass 
as a heat sink remains the same for the same blowing rate and, consequently, q w j/q w j 0 
decreases considerably with a small mixing coefficient. 

From the viewpoint of practical application, small coolant slot height seems to be 
preferable. The results of the present approach suggest that to reduce the increase of 
heating rate in the jet mixing region compared with that in the corresponding attached 
boundary-layer flow, a smaller value of Z is obviously preferable. Therefore, if the 
height of the coolant slot is fixed for the intended blowing rate, the stream wise location 
of the slot should be near the leading edge, where the thickness of the oncoming boundary 
layer is small. The reasons for this location can be summarized as follows: 

(1) The increase of mean heat -transfer rate for no coolant flow is decreased for 

smaller values of Z (fig. 13) 

(2) The coolant flow efficiency increases (fig. 16) 

(3) The length of the high heating region (length of separated region) decreases 

(fig. 11). 


Recompression Region 

The main difficulty which is encountered in numerical calculation for this region 
is evaluation of two parameters involved in the fundamental equations. These param- 
eters are (p and K, which in principle should be functions of k on the basis of the 
assumption of a one-parameter family for the boundary-layer state. Since there is 
neither analytical nor empirical information available for estimating these parameters 
at present, it is assumed that k = 0.03 (Crocco and Lees' value) and <p = 1.0 (flat- 
plate value). 

Figure 19 shows a typical example of the variation of local pressure, heat -transfer 
rate, and boundary- layer thickness in the recompression region without coolant mass 
injection. For this case the origin of the x-axis is taken at the exit of this region and 
6p denotes the boundary -layer thickness at the exit. The subscript °° indicates con- 
ditions of the external flow in the region of the fully developed turbulent boundary layer 
downstream of recompression, which are assumed to be the same as the oncoming flow 
conditions. As can be seen in figure 19, the region of interaction spreads in a width more 
than eight times the boundary- layer thickness in this region. In the case of shock-wave — 
turbulent-boundary-layer interaction, it is generally recognized that the shock wave dif- 
fuses to a width of the order of 10 times the boundary-layer thickness, and this fact seems 
to support the present result. 
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Figure 19 •- Local pressure, thickness of boundary layer, and heating 
rate in recompression region. =3-0; H = 0.3571; M = 0; 
N Pr,t = 1 -0* 


Thickness of the boundary layer becomes a minimum at a distance downstream of 
the entrance of this region and most of the recompression takes place before the minimum 
thickness is approached. Moreover, it is interesting to see in the figure that the maxi- 
mum heating rate occurs approximately at a point where the local pressure gradient 
becomes maximum. Furthermore, it should be remarked that this point does not coin- 
cide with the minimum boundary-layer thickness. These characteristics clearly indicate 
the physical fact that the local heat-transfer rate is more seriously influenced by pres- 
sure gradient (rate of change of external flow Mach number) than by local boundary-layer 
thickness. The decrease in local heating rate near the exit of the recompression region 
results from the combined effect of decrease in local pressure gradient and increase in 
boundary -layer thickness. 

The boundary -layer thicknesses at the entrance and exit of the recompression 
region are points of interest. Figure 20 shows the effect of coolant flow rate on the 
boundary-layer thickness at both edges of the recompression region. As can be seen in 
figure 20, boundary- layer thickness at the exit is rather small compared with that at the 
entrance for a very small coolant flow rate. This result is quite reasonable because 
recompression is strong enough for a low coolant flow rate to contract the boundary-layer 
thickness, even if a certain amount of mass may be transported into the viscous layer 
from the external stream, as has been already shown in figure 19. 
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(t>) Exit of recompression region. 

Figure 20.- Thickness of boundary layer. = 3*0; H c /H w = 1.0; 

Z = O.CA-2; N pr ^ t = 1.0. 

The boundary-layer thickness at the entrance of the recompression region increases 
almost linearly with the blowing rate and is rather insensitive to wall enthalpy. Although 
a lower wall enthalpy decreases the thickness of the viscous layer, it also elongates the 
separated region, which, in turn, increases the mass input into the viscous layer and 
increases its thickness at the exit of the jet mixing region. These two effects seem to 
cancel one another, and thus result in the insensitivity to the wall enthalpy of the boundary- 
layer thickness at the entrance of the recompression region. 

The difference between 6p and 6 r shown in figure 20 indicates an increment of 
boundary -layer thickness through the recompression region, but it is fairly small in spite 
of the mass transport through mixing. This condition may be due to a contraction effect 
through recompression. Figure 21 presents the effect of Z on boundary-layer thickness 
at both edges of the recompression region. Since Z increases either with an increase 
of oncoming boundary-layer thickness or with a decrease of slot height, the results shown 
in this figure are reasonable. 
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(a) Entrance of recompression region. 



(b) Exit of recompression region. 

Figure 21.- Thickness of boundary layer. 3^ = 3-0; 

Hw/H^ = 0.3571; H c /H w = 1.0; % r?t = 1.0. 

The width of shock-wave diffusion (width of recompression region) is another point 
of interest. Figure 22 shows the variation of width in shock-wave diffusion with the 
coolant flow rate, where the effect of wall enthalpy and the parameter Z are examined 
separately. In the case of shock-wave — laminar-boundary-layer interaction, it is gener- 
ally recognized that a lower wall enthalpy decreases the width of the shock-wave diffusion. 
The reason for this effect is that the extent of the subsonic region of the boundary layer 
near the wall, through which shock-wave disturbance may be felt upstream, decreases 
with decreasing wall temperature. The present results seem, at a glance, to be incon- 
sistent with the general characteristic just mentioned, since, as seen in figure 22, the 
width of the shock-wave diffusion increases with a decrease of the wall enthalpy even in 
the case of no-coolant flow. However, detailed examination of the velocity profile across 
the jet mixing layer reveals that the low velocity region of the velocity profile tends to 
diffuse far inside the dead-air region as the wall enthalpy decreases, and thus indicates 
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(b) = 3-0; Hw/Ho = 0.5571; Hc/H w = 1.0; N Pr?t = 1.0. 

Figure 22.- Width of shock wave diffusion. 

an increase in the subsonic region. The same might be true for the attached turbulent 
boundary layer although there is no rigorous proof for this statement. Thus, this char- 
acteristic seems to overcome the decrease of the subsonic region due to decreased wall 
enthalpy and results in a net increase of the subsonic region through which the shock- 
wave disturbance may propagate upstream. 

As the blowing rate increases, the effect of wall enthalpy on the width of the shock- 
wave diffusion seems, at a glance, to be pronounced. However, this condition is not due 
to the direct effect of wall enthalpy but mainly is due to differences in shock-wave strength 
for the same blowing rate. Since an increase in Z results in an increase of the 
boundary-layer thickness in the recompression region, the effect of Z on the width of 
the shock-wave diffusion is quite reasonable. (See fig. 22(b).) 

Figures 23 and 24 where the effects of wall enthalpy and Z are examined sepa- 
rately present the mean heat -transfer rate without coolant mass injection in the recom- 
pression region. In contrast to the result shown in figure 13, the mean heat-transfer 
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Figure 23.- Effect of enthalpy ratio on mean heat- transfer 
rate in recompression region, 1 ^,, = 3 .O; M = 0; 

Z = 0.042; N Pr t = 1.0. 
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rate in the recompression region decreases slightly with increased Z, as can be seen 
in figure 24. Clearly, the shock-wave strength decreases with increased Z, whereas 
the width of the shock-wave diffusion increases. This characteristic obviously indicates 
that the local pressure gradient decreases considerably with an increase in Z and, con- 
sequently, the maximum value of the local heating rate within the recompression region 
decreases to yield a net decrease of mean heating rate. 

Figure 25 in which the effect of wall enthalpy is examined shows the variation of the 
mean heat -transfer rate with the coolant flow rate. The mean heating rate decreases with 
an increase in coolant flow rate. However, since the mixing between the coolant flow and 
the oncoming flow has almost finished before the recompression region is approached, it 
seems that the coolant mass no longer acts as a heat sink in this region, but only con- 
tributes to increased boundary-layer thickness; therefore, the effectiveness of coolant 
flow, expressed as q w r 0 becomes higher compared with its value in the jet 
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Figure 25-- Variation of mean heat-transfer rate in re compress ion 
region with coolant flow rate for various enthalpy ratios. 

= 5-0; H c /H v =1.0; Z = 0.0^2; N pr ^ t = 1.0. 



Figure 26.- Variation of mean heat-transfer rate in recom- 
pression region with coolant flow rate. = 3-0; 

V H » = °- 3571; H =/ H w = 1 -°; N Pr,t = 1 -°- 

mixing region. Figure 26 shows the effect of Z on heating rate in the recompression 
region and indicates that the mean heat -transfer-rate ratio seems to be insensitive to Z 
in this region as well as in the jet mixing region. 

As has been mentioned previously, the local heat -transfer rate depends strongly 
upon the local rate of energy transport into the internal viscous layer from the external 
isentropic stream. This dependence can be interpreted to mean that the local heating 
rate is seriously influenced by the order of magnitude of the local mixing coefficient. In 
order to demonstrate this effect, the mean heat -transfer rate in the recompression region 
is plotted against the mixing coefficient in figure 27, where the shock-wave strength and 
wall enthalpy are fixed and no coolant mass injection is assumed. Figure 27 shows that 
the mean heat -transfer rate decreases considerably as the mixing coefficient decreases. 
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Figure 27-- Effect of mixing coefficient on mean 
heat-transfer rate in recompression region. 

M oo = 5.0; = 0.3571; Mx = 3.51^; 

M = 0; Z = 0.042; Np r? t = 1-0. 


In the preceding section, a qualitative discussion on heating rate of turbulent bound- 
ary layer has shown that 

q oc K6 -1 / 4 
H w 

The thickness of the boundary layer is proportional to the mass flux across the boundary 
layer, and the mass flux, in turn, is proportional to the mixing coefficient. Therefore, 
the relation may be reexpressed as 


Hence, 


q oc K' 

H w 


3/4 




(107) 


where q w ^ and q w g denote the mean heating rates corresponding to Kj and K^, 
respectively. Examination of the result shown in figure 27 clearly reveals that equa- 
tion (107) is a fair approximation. This result supports the validity of the discussion 
made in the preceding section on the characteristic behavior of the mean heat -transfer 
rate in the jet mixing region. 


Figure 28, where the effects of wall enthalpy and the parameter Z are examined, 
shows the variation of coolant flow Mach number at the slot exit at which the re compres- 
sion shock waves disappear. The critical coolant flow Mach number Me*, defined by 
the disappearance of the recompression shock wave, increases with either an increase in 
wall enthalpy or a decrease in Z. Reference 4 found experimentally that for isoenergetic 
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Figure 28.- Coolant flow Mach number at slot exit for 
disappearance of shock wave. = 3*0; H c /H w = 1.0; 

X = 1.0; N prjt = 1.0. 

injection (H c = H M ) and Z = 0.02, the disappearance of the recompression shock wave 
occurs before the choking condition of coolant flow is approached at the slot exit. This 
condition means that the critical Mach number remains subsonic. A rough evaluation 
shows that the critical coolant flow Mach number is about 0.65. This value is supported 
at least qualitatively by the experiment. 

Region of Fully Developed Turbulent Boundary Layer and Overall Results 

Reduced heat transfer in the region of a fully developed turbulent boundary layer 
results from the growth of boundary -layer thickness associated with the coolant mass 
injection upstream. Since the Mach number of the external stream is assumed to be kept 
constant at its oncoming value, a complicated flow phenomenon due to change in the exter- 
nal flow Mach number does not appear in this region. Moreover, a large amount of mass 
transport from the external stream to the internal dissipative flow in the upstream regions, 
which has been shown to have rather unfavorable effects on heating rate in the upstream 
regions, becomes a favorable factor in this region and results in reduction of local heat- 
transfer rate. This net reduction of overall heat transfer seems to be a characteristic of 
film cooling through a rearward-facing slpt. 
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Figures 29 and 30, where the effects of wall enthalpy and slot height are examined, 
respectively, show local effectiveness of film cooling in the region of the fully developed 
turbulent boundary layer. The local effectiveness of film cooling in this region decreases 
monotonically with increasing coolant flow rate. The reason may be understood through 
the subsequent argument. 


Since the local heat -transfer rate is approximately proportional to 6“*/^, 
easily shown that 


it is 


J w 

l W ,o 


r 


(108) 


where the subscript o denotes conditions of no coolant flow. As can be seen in fig- 
ures 20 and 21, the boundary-layer thickness Sp at the entrance of this region increases 
with increase of coolant flow, but the rate of increase becomes smaller as the blowing 
rate grows. These characteristics clearly indicate that the rate of increase of the 



Figure 29*- Local effectiveness of film cooling in region of fully 
developed turbulent boundary layer. = 3*0: Hr/H w = 1.0; 

z= 0.04 2; N RejD = icA ; x = 1.0; N prjt =i.o. ' 
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Figure 30.- Local effectiveness of film cooling in region 
of fully developed turbulent boundary layer. = 3-0; 

H w /lio = 0.3571; Hc/H w = 1.0; Z = 0.0^2; X = 1.0; 

N P r,t = 1-0. 

numerator in equation (100) is less than M except in the vicinity of zero coolant flow 
rate, where it is of the order of M, whereas the denominator is always of the order of 
M. Therefore, is finite at M = 0 and decreases monotonically as M grows. 

The effect of wall enthalpy on r ^ shown in figure 29 seems to be physically reason- 
able because the heat -transfer reduction is proportional to Hoo - H w . It is interesting to 
see in figure 30 that a smaller coolant slot height may be more effective. However, 
attention must be paid to the fact that the results shown in figure 30 do not indicate a pure 
effect of slot height, since the change of Np^ e p, the value of Z being kept unchanged, 
means a change in thickness of the oncoming boundary layer at the same time. It is clear 
that a decrease of Nr 6> p), Z being kept unchanged, means a decrease of 5 <*,, and indi- 
cates an upstream movement of the coolant slot location. This condition means a thin 
boundary layer at the entrance of the region of the fully developed turbulent boundary 
layer. In any event, the local heating rate decreases rather slowly, as is seen in the 
relation (108), with an increase of boundary- layer thickness due to coolant flow rate. 

It must be noted that in the region of the fully developed turbulent boundary layer, 
the magnitude of the local heat -transfer-rate reduction due to an increase of 
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boundary-layer thickness is not so large because q w °c However, since the region 

of integration used for obtaining r/f extends to infinity, the local effectiveness of film 
cooling in this region seems to become rather large. In fact, if an observation point for 
the local effectiveness of film cooling is taken at a finite distance downstream of the 
entrance of this region, then rjf will decrease considerably depending upon the location 
of the observation point. 

Figures 31(a) to 31(c) present the overall effectiveness of film cooling for wall 
enthalpy ratios of 0.5, 0.3571, and 0.1, respectively. The overall effectiveness of film 
cooling 77 t has an explicit physical meaning of total heat -transfer reduction pertinent to 
whole flow field downstream of the injection point. Therefore, it can be considered as a 
useful quantitative measure to indicate efficiency of the injected coolant mass. 



(a) NU = 3.0; Hw/H* = 0-5; H c /H w = 1.0; Z = 0.042; 
X = 1.0; N Pr>t = 1.0. 



M 

(t>) Moo - 3-0; Hy/Hoo = 0.3571; H c /H w = 1.0; Z= 0.042; 
X = 1.0; N p t = 1.0. 

Figure 31-- Overall effectiveness of film cooling. 
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(c) Mco =3-0; Hw/Hoo = 0.1; H c /h w = 1.0; Z = 0.042; X = 1.0; 

®Pr,t = I’®- 

Figure 31 •- Concluded. 

It is interesting to see in figure 31 that the magnitude of the overall effectiveness of 
film cooling becomes small compared with for low coolant flow rates, whereas it 
becomes larger than rjf as the blowing rate grows. This result clearly indicates the 
effect of the jet mixing and recompression regions. The mean heating rate in these high 
heating regions decreases monotonically with an increase in the coolant flow rate, but 
their length increases with blowing rate, for small blowing rates. These characteristics 
increase the net heat transfer to the wall within the high heating regions in the case of 
small coolant flow rate, and thus decrease On the other hand, as the blowing rate 

goes beyond a certain value, the mean heat -transfer rate within the high heating regions 
is still decreasing, whereas the total length of these regions becomes insensitive to cool- 
ant flow rate, so that the overall heating decreases. As a result, reaches its maxi- 

mum at a certain blowing rate which is dependent on wall enthalpy and N ReD . 

In order to confirm this point, the effect of wall enthalpy on r? T is examined by 
keeping other conditions constant and the result is shown in figure 32. The lower wall 
enthalpy seems to increase slightly the coolant flow rate at which Tj-p reaches a maxi- 
mum. Figure 33 presents an examination of the effect of the oncoming boundary-layer 
thickness on r) T . Since N Re D is kept constant, the increase of Z indicates an 
increase of 6*,. As is seen in figure 33, a small thickness of the oncoming boundary 
layer is favorable for increasing the overall effectiveness of film cooling, but the coolant 
flow rate for a maximum t? t is less sensitive to the oncoming boundary-layer thickness. 

Figure 34 shows the effect of the coolant slot height on the overall effectiveness of 
film cooling. The figure clearly indicates that a larger slot height results in a larger 
coolant flow rate at which r/j, becomes a maximum. Therefore, it may be concluded 
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Figure 3^*- Overall effectiveness of film cooling for 
selected conditions. = 3*0$ H w /Hoo = 0 . 3571 ; 

H c /H w - 1.0; X = 1.0; N pr ^ t = 1.0. 


that the blowing rate for maximum overall effectiveness of film cooling increases with 
increased slot height, and this trend becomes increasingly pronounced as the wall enthalpy 
decreases. 

It should be noted that there exists an optimum value of the coolant slot height for 
which the overall effectiveness becomes a maximum if the other conditions are fixed. In 
the regions of jet mixing and recompression, the coolant slot height determines the flow 
characteristics, mainly in a form of Z /= — ^ Y whereas in the region of fully 

\ D v^V 

developed turbulent boundary layer, it determines the flow characteristics in forms of Z 
and D itself. Since a complicated nonlinear coupling of these factors is an essential 
physical feature of the overall effectiveness of film cooling, it is not easy to give a com- 
prehensive argument on the physical reason for the existence of an optimum coolant slot 
height. For this reason, a simple discussion is given to assist in understanding the pres- 
ent result. 

Insofar as the present problem is concerned, the overall effectiveness of film 
cooling can be considered to be a function of Z, D, and M: 


t? t = f n (Z,D,M) 

Differentiation of 77-p with respect to D leads to 

d?7 T _ ( 5 n\ ( M n\ dZ 

dD \ 0 D/ Z + \ 8Z /D dD 


(109) 
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It is obvious from figure 31(b) that the first derivative on the right-hand side of equa- 
tion (109) is always negative. In the same way, the second derivative is always negative 
from the results shown in figure 33, whereas the third derivative is also negative. These 
characteristics suggest implicitly that might vanish at a certain finite value 

of D, and thus indicate that 77 , p has an extremum, which will be found to be a maximum 
from a rough examination. 


CONCLUDING REMARKS 

Heat transfer downstream of a two-dimensional rearward-facing step with coolant 
injection has been analyzed under the assumption that the oncoming boundary -layer thick- 
ness is small compared with the step height, and the coolant flow rate is so small that the 
velocity in the dead-air region is negligible. The whole flow field is divided into three 
regions: a jet mixing region, a recompression region, and a fully developed turbulent 

boundary layer far downstream. 

The flow field in the jet mixing region is shown to be determined uniquely for given 
conditions, and the pressure in this region increases with increased coolant flow. The 
length of the separated region also increases with increased blowing rate, is fairly sensi- 
tive to the oncoming flow conditions, and is rather insensitive to wall enthalpy. With no 
coolant flow, the mean heat-transfer rate within this region increases considerably com- 
pared with that in the corresponding attached boundary-layer flow. However, the mean 
heating rate decreases with the increased coolant flow. 

The increase in local heating rate compared with that in the corresponding flat- 
plate flow is shown to be inevitable in the recompression region, because of a strong 
adverse pressure gradient, and the local heating rate reaches a maximum at the point 
within this region where the local pressure gradient also is highest. The mean heat- 
transfer rate within this region decreases with increased blowing rate. However, in 
contrast to the jet mixing region, the decrease in mean heat-transfer rate is not large in 
this region. Thus, the jet mixing and recompression regions have a rather high heating 
rate compared with the corresponding flat-plate flow. 

The increased heating in these upstream regions is shown to be compensated suffi- 
ciently in the region of a fully developed turbulent boundary layer; as a result, the film 
cooling is fairly effective in reducing the overall heat transfer to the wall. It is shown 
that a smaller thickness of the oncoming boundary layer gives higher effectiveness of 
film cooling, and indicates that the upstream location of coolant slot is preferable. 

The overall effectiveness of film cooling, which increases with increased slot 
height, reaches a maximum at a very small coolant flow rate; this optimum blowing rate 
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depends slightly upon the wall enthalpy. The present approach further indicates that if 
the thickness of the oncoming boundary layer is fixed, there exists an optimum height 
of the coolant slot for which the overall cooling effectiveness is a maximum. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 29, 1970. 
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APPENDIX 


SOLUTION OF FUNDAMENTAL EQUATIONS IN JET MIXING REGION 

If u*/h* 2 = f(£) is assumed to be known, equation (9) can be written in a form 

+ (?+*’)—■ =0 (Al) 

d? 2 d? 

where the prime indicates differentiation with respect to £• Equation (Al) is an ordinary 
differential equation of the second order; thus, a general solution can be easily obtained as 

" - (A2) 


u* = C x f f d ? + C 2 
1 


where 


F(0 = exp 




(A3) 


and where Cj and C 2 are integration constants to be determined by the given boundary 
conditions. From equations (12) 


1 = C- 


if* 


1 0 I - + C 2 


0 = C lI 0 "f d? + C 2 


(A4) 


which yields 




C, = 


— *0 


d? 


(A5) 


J 


Therefore, the solution for u*(£) satisfying the given boundary conditions can be 
obtained as 


u*(0 = l yd? 


(A6) 
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APPENDIX 


The solution for h*(£) can be obtained similarly. Rewriting equation (10) yields 
the following ordinary differential equation of the second order, since u*(£) has been 
known from equation (A 6 ), 

f + ( Npr ’ t? * = " Npr ’ t(y ■ 1)M i 2f (tr) 2 (A7) 

Since 



a general solution to equation (A7) may be obtained as 


* 

h = C. 


I 


? F N Pr,t 


C F N Pr,t F 2 ' N Pr,t 


d? - /3 2 N p t (y - 1) M X 2 J E—1l J 

7 — oc I 


'-00 f J ~'"~ A J-oo f J 0 f 

where Cg and C 4 are integration constants. From equations (12) 

f » F Npr,t 2 2 F N Pr ,t f i F 2 - N Pr,t 

1 = C 3 J.„ — — dC - 13 Npr .‘ (r - 1)M l J.„ —T~ 1 f 


d ?7 d£ + C 4 (A 8 ) 


~\ 


d?7 d£ + C4 


> (A9) 


h c * = C 4 


which leads to 


Co = 


r 


p N Pr,t 


d? 


1 - h c * + ^ 2 N Pr t (r - l)M x 2 f 

7 ^ - OO 


00 F N Pr,t F 2-N Pr,t 
f Jn f 


d ?7 d| 


(A10) 


Thus, the solution for h (£) satisfying the given boundary conditions can be expressed as 


h*(?) = h c * - ^ 2 N Pr)t (y - 1)M X 2 J ^ F j T,t J ^ F f Pr,t dT? d£ + C 3 J 


2-N 


Pr,t 


? p^Pr ,t 


d? 

(All) 


where C 3 has been given by equation (A10). 
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